TO APPEAR IN THE ASTROPHYSICAL JOURNAL 
Preprint typeset using IATEX style emulateapj v. 5/2/11 


VALIDATION OF TWELVE SMALL KEPLER TRANSITING PLANETS IN THE HABITABLE ZONE 


GUILLERMO TonRES!?, Davip M. Kippinc!, FRANCOIS FRESSIN!, DOUGLAS A. CALDWELL?, JOSEPH D. TWICKEN?, SARAH 
BALLARDÓ?, NATALIE M. BATALHAÓ, STEPHEN T. Bryson®, DaviD R. CrARDI', CHRISTOPHER E. HENZE?, STEVE B. 
HowELL^, HOWARD T. Isaacson’, Jon M. Jenkins®, PHILIP S. MurmuEAD??, ELISABETH R. NEWTON!, ERIK A. 
PETIGURA?, THOMAS BARCLAY’, WILLIAM J. Borucki®, Justin R. Creppe!!, MARK E. EvERETT, ELLIOTT P. HoncH??, 
ANDREW W. Howarp", REA KOLBL?, GEOFFREY W. Marcy’, SEAN MCCAULIFFP, AND ELISA V. QUINTANA? 


To appear in The Astrophysical Journal 


ABSTRACT 


We present an investigation of twelve candidate transiting planets from Kepler with orbital periods 
ranging from 34 to 207 days, selected from initial indications that they are small and potentially in 
the habitable zone (HZ) of their parent stars. Few of these objects are known. The expected Doppler 
signals are too small to confirm them by demonstrating that their masses are in the planetary regime. 
Here we verify their planetary nature by validating them statistically using the BLENDER technique, 
which simulates large numbers of false positives and compares the resulting light curves with the 
Kepler photometry. This analysis was supplemented with new follow-up observations (high-resolution 
optical and near-infrared spectroscopy, adaptive optics imaging, and speckle interferometry), as well as 
an analysis of the flux centroids. For eleven of them (KOI-0571.05, 1422.04, 1422.05, 2529.02, 3255.01, 
3284.01, 4005.01, 4087.01, 4622.01, 4742.01, and 4745.01) we show that the likelihood they are true 
planets is far greater than that of a false positive, to a confidence level of 99.73% (3e) or higher. 
For KOI-4427.01 the confidence level is about 99.2% (2.60). With our accurate characterization of 
the GKM host stars, the derived planetary radii range from 1.1 to 2.7 Rẹ. All twelve objects are 
confirmed to be in the HZ, and nine are small enough to be rocky. Excluding three of them that have 
been previously validated by others, our study doubles the number of known rocky planets in the HZ. 
KO1-3284.01 (Kepler-438 b) and KOI-4742.01 (Kepler-442 b) are the planets most similar to the Earth 
discovered to date when considering their size and incident flux jointly. 

Subject headings: methods: statistical — planetary systems — stars: individual (KOI-0571 — Kepler- 
186, KOI-1422 = Kepler-296, KOI-2529 = Kepler-436, KOI-3255 = Kepler-437, 
KOI-3284 — Kepler-438, KOI-4005 — Kepler-439, KOI-4087 — Kepler-440, KOI- 
4427, KOI-4622 = Kepler-441, KOI-4742 = Kepler-442, KOI-4745 = Kepler-443) 
— techniques: photometric 
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1. INTRODUCTION 


Over the duration of its four-year mission the Kepler 
spacecraft has enabled the identification of several thou- 


sand candidate transiting planets ( Kepler Objects of In- 
terest, or KOIs; 2011} [Batalha et al. 
, and many more continue to be found 


from reanalysis of the original data with increasingly so- 
phisticated methods. Only a tiny fraction of these can- 
didates have been “confirmed” in the traditional sense 
of having had their masses measured, either spectro- 
scopically or by modeling their transit timing variations 
(TTVs). Hundreds of others, mostly in multiple systems 
(*multis"), have recently been shown statistically to have 
a very high chance of being true planets 
bord [Rose et aL P014). even if their masses are not cur- 
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rently known. 

Public curiosity and scientific interest have motivated 
efforts in the last few years to find and confirm rocky 
planets similar in size to the Earth that are orbiting in 
the so-called habitable zone (HZ) of their parent stars, 
usually taken in this context to be the region in which 
water on the surface can be in a liquid state. Recent es- 
timates by [Petigura et al] suggest the rate of oc- 
currence of Earth-size planets in the HZ of Sun-like stars 
(which they defined broadly as corresponding to an inci- 
dent flux between 4 and 4 times that of the Earth) may 


be as high as 22+8% (for other estimates see 
Foreman-Mackey et al] (2014), 


but very few of these have actually been confirmed. 
These two conditions for habitability — rocky nature 
and suitable location — are sometimes more challenging 
to establish unambiguously than it may seem, as they 
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require an accurate knowledge of the host star prop- 
erties, something that is not always trivial to achieve 
for the faint Kepler targets. 


Definitions of the HZ 


radius at which planets transition from being rocky to 
being dominated by thick hydrogen/helium envelopes, 
which presumably makes them unsuitable for life as 
we know it. The interface is currently thought to 
be between 1.5 and 2 Ha (e.g. 


and it is possible that some planets that are slightly 
larger may still be rocky. 

These difficulties notwithstanding, a handful of tran- 
siting planets that appear to meet these conditions have 
already emerged from the Kepler sample. Among those 
considered at the time of publication to be smaller than 
2.5 Ra and in or near the HZ, the first to be announced 
was Kepler-22b (Borucki et a1][2012), with a radius of 
Rp = 2.38 + 0.13 Ra. Others followed including Kepler- 
6lb (2.15 + 0.13 Ro; 2013), Kepler-62e 
and f (1.61 + 0.05 Rg and 1.41 € 0.07 Ra; 
(2013), Kepler-69c (1.717934 Ra; [2013), 
and recently Kepler-186 f (1.1140.14 Rẹ ;[Quintana et al] 
[2014). All of these were “validated” statistically using 
techniques such as BLENDER 
by examining the likelihood of dif- 
ferent false positive scenarios to measure the confidence 
in the planetary interpretation. Some of them are in 
multi-planet systems; as mentioned above,|Lissauer et al. 
have developed a statistical argument by which 
they demonstrate that most candidates in multis are 
very likely to be true planets. Based on this reason- 


ing and a careful examination of follow-up observations, 
reported an additional small planet 
(Kepler-296 f, 2.31 + 0.12 He ) that appears to be in the 
HZ, according to (2014). Taking advantage 
of the same statistical framework for validation, the lat- 
ter authors announced one more small HZ planet in 
a multi-planet system: Kepler-174d (2.19 + 0.13 Ra). 
Other small HZ candidates have been proposed (see, e.g., 
2012: 
Gaidos |2013; [Mann et al. 
but have yet to be investigated in detail 
and confirmed. A few additional examples also believed 
to be in this special group were later shown to fail one or 
both requirements (size and location), either because of 
incorrect stellar parameters (originally based on photom- 
etry and later improved with spectroscopy) or because of 
the discovery of other stars in the photometric aperture 
that dilute the transits and change the inferred planetary 
radius (e.g., [Star et al][2014]. 

In this paper we investigate a sample of twelve Kepler 
candidates identified on the basis of preliminary analy- 
ses indicating that they are small and likely to be in or 
near the HZ of their parent stars. Our goal is to val- 
idate them as bona-fide planets, and to confirm their 
key properties pertaining to habitability. We describe 
our efforts over the past year to obtain the follow-up 
observations necessary to robustly characterize the host 


stars and to validate the signals. In the interim, three of 
them that are in systems with five candidates each (KOI- 
0571.05 — Kepler-186 f, KOI-1422.04 — Kepler-296 f, and 
KOIL1422.05 = Kepler-296e) were the subject of re- 
cent studies by others. KOI-1422.04 and KOI-1422.05 
were announced as planets by [Lissauer et al] and 
using their statistical framework for 
validation of multis “in bulk”, though only the first was 
mentioned as being in the HZ. KOI-0571.05 was validated 
independenti byl@uumtaas ct E. These acthors 
relied in part on the fact that false positives are much less 
common in multis, and on the validation of the other four 
planets in the same system by the previous authors, in or- 
der to reach a sufficiently high level of confidence for the 
fifth candidate. Most of the other targets in our sample 
are considerably more difficult to validate because they 
tend to have long orbital periods and do not belong to 
multiple systems. Nevertheless, as we show below, the 
application of the BLENDER technique used in many of the 
previous discoveries is able to achieve validations here as 
well. Our work has now essentially doubled the number 
of known small HZ planets. 

Our paper is organized as follows: Sect. [2] and Sect. 
describe our sample of candidates and the Kepler pho- 
tometry we use. The follow-up observations including 
high-resolution imaging, optical spectroscopy, and near- 
infrared spectroscopy are presented in Sect. |4] where we 
also report the complicating discovery of nearby stellar 
companions to four of the targets. These close compan- 
ions not only dilute the transits but also bring ambiguity 
as to the precise location of the planets in these sys- 
tems. A description of the analysis of the flux centroids 
that addresses some of those companions is given here as 
well. Sect. [5] follows with a determination of the stellar 
properties. The formal validation of our candidates is 
the subject of Sect. [6] after which we proceed with the 
transit light curve fits to derive the planetary parameters 
(Sect. [7). In Sect. [8] we apply the powerful technique of 
asterodensity profiling to investigate whether the planets 
that are in multiple systems orbit the same star, to ex- 
tract eccentricity information from the light curves of our 
targets, and for the four host stars with close stellar com- 
panions, to address the possibility that the planets orbit 
the companions rather than the intended targets. The 
topic of habitability is examined in Sect. D} We conclude 
with a discussion of our findings and final remarks. 


2. CANDIDATE IDENTIFICATION 


Candidate transiting planets of special interest for 
this study were initially identified as ‘threshold crossing 
events’ (TCEs; objects with a 7.1c or greater multiple- 
event transit detection statistic; see Dorkins et allo 
from a transit search conducted using nearly three 
years of Kepler data gathered from quarters 1 to 12, with 
version 8.3 of the Science Operations Center (SOC) pro- 
cessing pipeline {Tenenbaum et all D013). The analysis 
of three years of data gave us the first chance to de- 
tect three transits from planets in habitable-zone orbits 
around stars like the Sun, as well as to significantly in- 
crease the signal-to-noise ratio (SNR) for transit signals 
from habitable-zone orbits around cooler stars. In order 
to investigate the smaller planets potentially in the HZ, 


we selected for further scrutiny the TCEs that met the 
criteria Rp < 2 Ra and Teq < 303K, as well as a transit 
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SNR > 7. Here Teq represents the planetary equilibrium 
temperature, calculated assuming an Earth-like albedo, 
and the stellar properties relied on information from the 
Kepler Input Catalog (KIC; [Brown et_al2011)). We used 
data products from versions 8.3 and 9.0 of the SOC Data 
Validation (DV) pipeline module for the initial triage and 
vetting (for an overview of this process, see Wir et al] 
[20109]. We compared the results of our vetting with the 
scores from early runs of an automated vetting code be- 
ag developed at the iras Mecake albo Ud. Fred 
the 18,407 TCEs that passed through DV, we identified 
385 TCEs around 338 unique stars that met our selec- 
tion criteria. From this list of 385 TCEs we identified ten 
that passed flux and centroid vetting criteria and that 
had low false-alarm probabilities from the auto-vetter. 
During this vetting process, the SOC completed a tran- 
sit search of Kepler data from quarters 1-16 using ver- 
sion 9.1 of the processing pipeline. With the addition of 
three quarters of data and the improvements in the ver- 
sion 9.1 DV products, we promoted an additional TCE 
(KOI-4742.01), which had marginal signal-to-noise ratio 
in the 12-quarter run, and identified KOI-1422.05, the 
fifth KOI and second habitable-zone candidate around 
this star. The complete list of targets selected for this 
study is given in Table [] 

We note that there is significant contamination of the 
TCE list for periods near the orbital period of the Ke- 
pler spacecraft (372 days) caused by instrumental im- 
age artifacts on specific detector channels, as has been 


noted by (2013). The contamination 
is due to moiré pattern noise injected during readout 
of some of the detector chains 
[Kolodziejczak et al] 2010). The moiré noise generates 
a host of single-event transit-like signals for targets on 
the noisy channel. When folded at a period near one Ke- 
pler year these signals can give multiple-event statistics 
above the 7.1o detection threshold, as the same target 
star falls on the same noisy channel. We have not in- 
cluded any KOIs from the moiré pattern noise channels 
in our analysis here. Two of the KOIs (KOI-2529.02 and 
KOI-4005.01) do show some of their transits on detec- 
tor channel 56, which has increased noise relative to the 
other channels on which these targets fall. However, we 
included these KOIs in our sample because both show 
transits on other channels and both passed all of our 
vetting checks. 


3. KEPLER PHOTOMETRY 


For the analysis in this paper we have made use of 
the publicly available Kepler data for each KOI from the 
(MAST)IP] The 
observations span Kepler quarters 1-17, corresponding to 
a period of four years beginning in 2009 May. The down- 
loaded data were made available as part of Data Release 
23 and were processed using SOC Pipeline version 9.1. 
While long-cadence (29.4 min) data were available for all 
KOIs, short-cadence (58.9s) data were only available for 
KOI-0571 and KOI-1422. In the following we describe 
our further processing of the simple aperture photome- 
try (SAP) time series used for the detailed light-curve 


16 All DV results are available from the NASA Exoplanet 


Archive: http: //exoplanetarchive.ipac.caltech.edu|. 


17 https://archive.stsci.edu/index.html . 


analysis of all KOIs. Later in Sect. [5.1] and Sect. [6-1] we 
describe the slightly different processing of the photom- 
etry used to extract stellar rotation information and for 
the BLENDER analysis. 

Prior to making use of the SAP measurements we re- 
moved instrumental and stellar photometric variability 
in the Kepler data that can disturb the transit light 
curve profile. This process was conducted in two steps: 
pre-detrending cleaning, and long-term detrending. The 
pre-detrending cleaning was carried out independently 
for each transit event, restricted to plus or minus half 
an orbital period surrounding the mid-transit time. We 
visually inspected each epoch and removed any charge- 
trapping ramps, flare-like events, and instrumental dis- 
continuities in the data. We made no attempt to cor- 
rect these artifacts and simply excluded them from the 
photometry manually. We then removed all of the tran- 
sit signals of the other known candidates in each system 
(for stars having multiple candidates) within +0.6 transit 
durations centered on the mid-times, using the reported 
ephemerides. Finally, we cleaned the data of 30 outliers 
from a moving median smoothing curve with a 20-point 
window. 

Next we removed long-term trends, which can be due 
to instrumental effects such as focus drift, or stellar ef- 
fects such as rotational modulations. For this task we 
used the Cosine Filtering with Autocorrelation Mini- 
mization (CoFiAM) algorithm, which was specifically de- 
veloped to protect the shape of a transit light curve. 
We briefly highlight the main features of CoFiAM and 
the specific inputs used for this analysis, and we di- 
rect the reader to (2013) for a more de- 
tailed description of the technique. It is essentially a 
Fourier-based method that removes periodicities occur- 
ring at timescales greater than a pre-designated “pro- 
tected timescale". In this work, we selected three times 
the known transit duration to be the protected timescale, 
which ensures that the transit profile is not distorted 
in frequency space. CoFiAM does not directly attempt 
to remove high frequency noise, since the Fourier trans- 
form of a trapezoidal-like light curve contains significant 
high frequency power (Waman 2019). It is able to 
explore many different harmonics by trying longer pro- 
tected timescales than the nominal choice (we capped 
the maximum number of harmonics at 30) and evalu- 
ate the autocorrelation at a pre-selected timescale (we 
used 30 minutes) locally surrounding each transit. From 
these different harmonics, we selected the harmonic order 
that minimizes this local autocorrelation, as quantified 
using the Durbin-Watson statistic. This *Autocorrela- 
tion Minimization" component of CoFiAM provides opti- 
mized data for subsequent analysis. For each KOI we 
defined the local transit data as being within six tran- 
sit durations either side of the mid-transit time, in order 
to provide an ample out-of-transit baseline. These local 
data were divided through by the final CoFiAM function 
and then stitched together to form our final light curve 
for analysis. 


4. FOLLOW-UP OBSERVATIONS AND CENTROID 
ANALYSIS 


4.1. High-resolution imaging 


'The presence of other stars blended with our targets 
is a potential concern, as these companions could be or- 
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Table 


1 


Sample of KOIs in this study. 


Candidate Name KID 
KOI-0571.05  Kepler-186 f 8120608 
KOI-1422.04 Kepler-296f 11497958 
KOI-1422.05 Kepler-296e 11497958 
KOI-2529.02  Kepler-436 b 8463346 
KOI-3255.01  Kepler-437 b 8183288 
KOI-3284.01  Kepler-438b 6497146 
KOI-4005.01  Kepler-439 b 8142787 
KOI-4087.01 Kepler-440 b 6106282 
KOI-4427.01 4172805 
KOI-4622.01 Kepler-441b 11284772 
KOI-4742.01 Kepler-442 b 4138008 
KOI-4745.01 Kepler-443b 11757451 


Note. — Columns after the first indicate the Kepler planet designation, 


Kp b 
(mag) (deg) 
14.005 82 
15.921 18.0 
15.991 180 
15.856 . 6.5 
14.852 9.0 
14.467 18.6 
14.560 19.9 
15.134 15.8 
15.645 82 
15.142 19.1 
14.976 15.0 
15.891 17.0 


Period 
(days) 


129.94 
63.34 
34.14 
64.00 
66.65 
35.23 

178.14 

101.11 

147.66 

207.25 

112.31 

177.67 


Depth 
(ppm) 


540 
1070 
850 
1150 
520 
400 
700 
1010 
1220 
960 
560 
930 


Kepler identification number, brightness in the Kepler passband, Galactic 
latitude, orbital period, and transit depth in parts per million relative to the 
out-of-transit stellar flux. For consistency in this paper we will refer to all 
objects by their original KOI names throughout, even though the validations 
described later earn them the official Kepler planet designations listed in the 


second column. 


bited by another object causing eclipses that may be the 
source of the transit signals we observe. Even if they are 
not eclipsing, the light of the companion stars will atten- 
uate the signal and lead to a bias in the planetary radius 
determination if the effect is not accounted for. 

Images from the J-band UK Infrared Telescope sur- 
vey (UKIRT; [Lawrence et alJ[2007]) have indeed revealed 
companions to three of our targets (KOI-2529, 3284, and 
4427) that are close enough to fall within the photo- 
metric aperture of Kepler (i.e., within a few arc sec- 
onds), although as we describe later they are ruled out 
as the source of the transits by our centroid analysis in 
Sect. [4.2] We list these companions in Table] based on 
information taken from the Kepler Community Follow- 
up Observing Program (CFOP) Web sit} Additional 
companions to three other KOIs are seen in the UKIRT 
images at larger separations (and are also excluded by 
the centroid analysis), but these stars are all very faint 
(Am = 5), and have a negligible impact on the inferred 
planetary sizes. 

The UKIRT images have a typical seeing-limited res- 
olution of about 0"8 or 0"9. To explore the inner re- 
gions around our targets beyond the reach of UKIRT we 
observed them with near-infrared adaptive optics ( (AO) 
in the J (1.248 m and K' 2.12 filters using a 
NIRC2 i aee inima (Wizinowich ct aiaa ican 
on the KeckII, 10m Temm KOI-0571 and 
KOI- 3255 were observed in August 2012, and the rest 
in August 2013 as part of a general infrared AO sur- 
vey of KOls (e.g. Adams ct all 201% [Rowe et al] 2014 
ERES FT . KOI-1422 was not observed with AO 
for this project as it has been the target of a separates ef- 

Gilliland et al] 2013 Star et al. 


fort using HST (Gilliland et al. Star et al. 


those authors found it to have a close MUTET see 
below). For KOI-0571, 3255, 2529, 3284, and 4742 the 
targets themselves were used as natural guide stars; for 
KOI-4005, 4087, 4427, 4622, and 4745 we used the laser 
guide-star AO system. In all cases the observations were 
obtained in a 3-point dither pattern to avoid the lower 
left quadrant of the NIRC2 array, which displays elevated 


18 nttps://cfop.ipac.caltech.edu/home/ . 


Table 2 
Close companions to target stars. 


Ang. sep. P.A. Mag. diff. 

Star (^) (deg) (mag) Source 
KOI-1422 0.220 216.3 AR=1.72> Speckle 
KOI- Kp 0.402 350.4 AK'—4.69 Keck AO 

5.05 42.6 AJ-—3.75 UKIRT 
KOT- 3255 0.180° 336.5 AK’=0.05 Keck AO 
KOI-3284 0.438* 193.2 AK'—2.03 Keck AO 
ee 4.109 359.7 AJ —2.90  UKIRT 
KOI-4427 4.76 274.8 AJ = 2.87 UKIRT 


a For targets with both a close and a wide companion the 
magnitude differences for the wide companions are relative 
to the total brightness of the inner pair. 

b Also observed at 880nm (approximately Sloan z band), 
giving Az = 1.62. The separation and P.A. in the table are 
the average of the two speckle bands. Similar results were 


reported by [Star et al] (2014) from HST observations. 

© This companion was also detected in our speckle imaging 
observations. 

4 This companion was also detected in the U BV survey 


of[Everett et al] (2012), with magnitude differences AB = 
1.80 and AV = 2.01. The separation and P.A. in the table 
are the average of the two surveys. 


noise. Five images were collected per dither pattern po- 
sition, each shifted 0'5 from the previous dither position 
to enable the use of the source frames for the creation of 
the sky image. 

In general the NIRC2 imaging detects all of the sources 
found in the UKIRT J-band imaging within 5" of each 
target. The NIRC2 array has 1024 x 1024 pixels with 
a scale of about 10maspix-!, and a field of view of 
1071 x 1071. Each frame was dark-subtracted and flat- 
fielded, and the sky frames were constructed for each 
target from the target frames themselves by median fil- 
tering and co-adding the 15 or 25 dithered frames. In- 
dividual exposure times varied depending on the bright- 
ness of the target, but were long enough to obtain at 
least 5000 counts per frame (NIRC2 has a gain of 4 elec- 
trons per DN); frame times were typically 1-30 seconds. 
Data reduction was performed with a custom set of IDL 
routines. Close companions were detected around KOI- 
2529, 3255, and 3284 (see Figure[I], two of which also 
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Figure 1. Images in the K’ band from our Keck AO imaging 
observations showing close companions found around KOI-2529, 
KOI-3255, and KOI-3284 (see also Table[2). 
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Figure 2. Sensitivity curves from our Keck AO imaging observa- 
tions in the (2.124 um) K’ band. Curves correspond to the KOIs 
as labeled, top to bottom. 


have wider companions seen in the UKIRT images. We 
report the relative positions and brightness of these close 
companions also in Table [2] 

Point source detection limits were estimated in a se- 
ries of concentric annuli drawn around each target star. 
'The separation and widths of the annuli were set to the 
full width at half maximum (FWHM) of the target point 
spread function (PSF). The standard deviation of the 
background counts was calculated for each annulus, and 
the 5c limits were determined for each annular region 
(see also [Adams et al]2012). The PSF widths for our 
Keck images were typically found to be 4—5 pixels, cor- 
responding to 0704-0705 FWHM. Typical contrast levels 
from the images are 2-3 magnitudes at a separation of 1 


FWHM, 7-8 magnitudes at 54- FWHM, and deeper past 
10 FWHM. Sensitivity curves in the K’ band for each 
of the targets are shown in Figure[2] and extend to an- 
gular separations of approximately 3", set by the dither 
pattern overlap. We supplemented these measurements 
with similar sensitivity estimates made from the UKIRT 
images available on the CFOP Web site, which reach to 
much wider separations but are not as deep. 

Additional high-resolution imaging observations for 
four of our targets (KOI-0571, 1422, 3255, and 3284) 
were obtained with the Differential Speckle Survey In- 
strument (DSST; [Horch et al] (2009; 
mounted on the 8.1m Gemini-N telescope during the 
time period UT 2012 July 25-28. The seeing conditions 
were generally 075-0"8 during those nights. The large 
aperture of the Gemini-N primary mirror allows for high 
SNR observations at faint magnitudes and achieves a 
very small diffraction limit for visible-light observations. 
'The speckle camera provides simultaneous observations 
in two filters. For all Gemini-N observations discussed 
here we used filters centered at 692 nm and 880 nm (cor- 
responding very roughly to the Cousins R and Sloan z 
bands), with band widths of 40nm and 50nm, respec- 
tively. Although full details of the observing protocols 
have been given previously by Horch et a1] (2012), we 
summarize them here for completeness. 

For bright targets observed with this instrument (V « 
12) a single sequence of 1000 speckle frames of 60 ms du- 
ration each has usually been found to be sufficient, for a 
total of 3-4 minutes of observing time. The much fainter 
sources discussed here (15 « V « 17) required up to 
eight, 1000-frame sequences that were later combined, or 
about 30 minutes of on-source exposure time each. Cal- 
ibration images and sequences were taken as is usual to 
allow us to measure and define dispersion effects, position 
angles, faint limit sensitivity, photometric performance, 
and for point source reconstruction. 

The speckle observations revealed no stars in the vicin- 
ity of KOI-0571, but KOI-1422 was found to have a 
close companion (seen in both filters) that was also 
subsequently detected in the HST observations by 
(2015). Its brightness and position rel- 
ative to the primary are given in Table[2] Close compan- 
ions were also found around KOI-3255 and KOI-3284, 
which are the same as detected in our AO imaging with 
Keck. 

For estimating our detection threshold for faint com- 
panions we considered each local maximum in the re- 
constructed image as a potential stellar source, and de- 
termined the statistics of these peaks (i.e., their average 
values and standard deviation) as a function of angu- 
lar distance from the primary star. We then adopted a 
conservative 5e threshold for the detection of any com- 
panion stars. The details of the calibration procedures 
mentioned above and the precision that may be obtained 
are described by [2012). The sensitiv- 
ity curves for the four Kepler targets observed here with 
speckle imaging are shown in Figure [B] 


4.2. Centroid analysis 


One method to identify possible false positives due to 
background eclipsing binaries measures the location of 
the transit signal relative to each KOI host star via dif- 
ference imaging. For each quarter, a difference image 
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Figure 3. Similar to Figure[2]for the four targets observed with 
speckle interferometry at the Gemini North telescope. Curves are 
shown for the (692 nm) R-band observations, and correspond to 
the KOIs as labeled, top to bottom. 


may be formed by subtracting an average of in-transit 
pixel values from out-of-transit pixel values. If the tran- 
sit signal is due to a stellar source, the difference image 
will show that stellar source, whose location is deter- 
mined by Pixel Response Function centroiding (PRF; see 
[Bryson et abdo) This location may then be compared 
to the location of the host star, determined from a PRF- 
fit centroid of an average out-of-transit image (assuming 
the host is well isolated). Repeating this for each quar- 
ter in which transits occur gives a collection of offsets 
of the transit source from the host star. The average of 
these quarterly offsets defines the final offset of the tran- 
sit source from the host star. This metric is reported for 
each KOI on the NASA Exoplanet Archive. If the offset 
is greater than three times its uncertainty, there is cause 
for concern that the signal may be coming from a back- 
ground source rather than the KOI itself. In that case we 
would conclude that the likelihood that the transit signal 
is due to a planet around the KOI host star is low (for de- 
tails, see[Bryson et al|2013). Any transit source location 
less than 3c from the host star is considered statistically 
indistinguishable from the host star. This 3e exclusion 
radius depends on the KOI rather than the host star, as 
it is primarily a function of the signal-to-noise ratio of 
the KOF's transit signal. 

Table[3]lists the average offset of the transit source for 
each KOI from the KOI host star, as reported on the Exo- 
planet Archive, as well as the 3c exclusion radius for each 
KOI. All KOIs have offsets that are less than three times 
their uncertainties, with the exception of KOI-4427.01. 
Therefore, on centroid grounds alone KOI-4427.01 has a 
non-negligible chance of being a false positive. 

The companions at 4” to 5" from KOI-2529, 3284, and 
4427 listed in Table [2] are significantly farther than 3c 
from the respective observed transit positions, so they 
are ruled out as possible sources of the transit, as are 
any wider companions detected in the UKIRT images. 
On the other hand, the close (< 0/5) companions to 
KOI-1422, 2529, 3255, and 3284 cannot be ruled out by 
centroid analysis alone. We address this issue later. Ad- 
ditional closer (unresolved) companions may of course 
also be present. 


Table 3 
Centroid results for the candidates. 


Offset from Offset from 3c exclusion 
Candidate star (”) star (c) radius (") 
KOI-0571.05 0.74 + 0.82 0.90 2.45 
KOI-1422.04 0.34 + 0.34 1.00 1.03 
KOI-1422.05 0.54 + 0.33 1.67 0.98 
KOI-2529.02 0.08 + 0.64 0.13 1.92 
KOI-3255.01 0.55 + 0.19? 2.82 0.58? 
KOI-3284.01 0.60 + 0.94? 0.64 2.82? 
KOI-4005.01 0.23 + 0.61 0.38 1.84 
KOI-4087.01 0.11 + 0.41 0.26 1.22 
KOI-4427.0n ^ 0.83 + 0.24» 3.47 0.72 
KOI-4622.01 0.92 + 0.52 1.76 1.57 
KOI-4742.01 1.19 + 1.17 1.01 3.52 
KOI-4745.01 0.76 + 0.36 2.09 1.08 


? The presence of relatively bright close companions to KOI- 
3255 and KOI-3284 can cause biases in the centroid measure- 
ments. These biases are expected to be smaller than the sep- 
aration of the close companion, but we have nevertheless con- 
servatively doubled the uncertainty found on the Exoplanet 
Archive. 

b The Exoplanet Archive reports the offset of KOI-4427.01 as 
0"94-- 1" 57, but there is a single large outlier due to a quarter 
with a noisy difference image. The values reported here are 
after removing the outlier. 


The centroid offsets in Table [B| on the basis of which 
we rule out the 4—5" companions, were computed on the 
assumption that the host stars do not have closer (< 4”) 
companions of comparable brightness. However, as de- 
scribed in Sect. [1|] this is not the case for KOI-3255 
and KOI-3284. The relatively bright close companions 
for these two KOIs will introduce error in the centroid 
measurement of the KOI star position. (Such error does 
not arise for the difference-image centroids measuring the 
transit source position because, assuming negligible vari- 
ability, the difference image shows only the star hosting 
the transit signal. To account for this unknown error, 
we have chosen to conservatively double the offset un- 
certainties for KOI-3255 and KOI-3284, although we do 
not expect the error to be larger than the separation of 
the companion. The uncertainties in the table already 
include this conservative doubling. The companion for 
KOI-2529 is sufficiently dim that the centroid error due 
to this companion should be negligible. 


4.3. Optical spectroscopy 


In order to characterize the target stars and search for 
additional companions, we acquired high-resolution spec- 
tra for the eleven stars in our sample with the Keck I 
Telescope on Mauna Kea (HI) and the HIRES spectrom- 
eter (Vout et alll1994) using the standard setup of the 
California Planet Search (CPS; see [Howard et allDO10). 
The wavelength coverage is approximately 3800-8000 A, 
recorded on three contiguous CCD detectors. Obser- 
vations were made over the course of four nights (UT 
2012 July 26 and September 22, and 2013 June 7 and 
28) under clear skies, with a typical seeing of ~1”. Ex- 
posure times of about 15 minutes resulted in signal-to- 
noise ratios of 12-25 per pixel at a mean wavelength 
near the Mg Ib triplet (5150-5200 A). The exception was 
KOI-3284, which was observed for 45 minutes yielding 
a SNR of 45 per pixel. All spectra were taken using 
the C2 decker (0786 x 14"), giving a resolving power of 
R ~ 60,000. The 14" long slit allows for subtraction of 
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the sky background that is superimposed on the stellar 
spectrum. Sky subtraction was performed by measuring 
the sky background light on a wavelength-by-wavelength 
basis. 

For the reduction of the raw spectra we used the stan- 
dard CPS pipeline. Images from each of the three HIRES 
CCDs were independently reduced and subjected to sky 
subtraction, flat fielding, and cosmic ray removal. The 
pixel columns at each wavelength were then summed, 
providing photon counts as a function of wavelength for 
each pixel. Consistent wavelength solutions were ensured 
by aligning a carefully chosen set of Thorium-Argon 
emission lines onto the same pixels at the beginning of 
each night’s observations. 

The radial velocity (RV) of each star was measured 
using the A-band and B-band telluric line features as 
wavelength fiducials. We measured the placement of the 
stellar absorption lines relative to these features, and 
then referenced them to stars of known radial velocity 
(Chubak et a1][2012). The final radial velocity measure- 
ments in the frame of the solar system barycenter are 
accurate to +0.1kms~!, and are given in Table 4] with 
other details of the observations. 

To aid in eliminating some of the false positives that 
might be causing the transit-like signals in our targets, 
we searched each of our spectra for evidence of a sec- 
ond set of lines from another star possibly falling on 
the spectrometer slit during our observations (i.e., at 
angular separations smaller than about 0743, or half 
of the slit width). We first cross-correlated each spec- 
trum against a large library of stellar spectra obtained 
with the same instrument and setup (see Sect. [5) in or- 
der to identify the best match. After flux normaliza- 
tion and placement of the target star and library stars 
onto a common wavelength solution, we subtracted this 
best match from the observed spectra, and the residu- 
als were subsequently cross-correlated against the same 
set of library stars (Kolbl et al]|2015). A resulting peak 
in the cross-correlation function would be evidence of a 
second star contaminating the spectrum. For our targets 
this method is sensitive to physically associated as well 
as unrelated (background/foreground) companion stars 
so long as the relative RV of the two objects is greater 
than about lOkms-!. For smaller velocity differences 
the technique loses sensitivity due to imperfect subtrac- 
tion of the primary star. 

We detected no companions to any of the 11 stars 
in our sample down to about 196 of the brightness of 
the primary (corresponding to a magnitude difference 
A Kp = 5) for targets brighter than Kp = 15, and down 
to 2% (AKp e 4.2) for our fainter targets. This indicates 
that at least two of the close companions found in our 
Keck AO imaging that are bright enough to have been 
detected spectroscopically, those around KOI-3255 and 
KOI-3284, must have radial velocities within 10kms-! 
of the main star, which strongly suggests they are phys- 
ically associated with the target. 


4.4. Near-infrared spectroscopy 


For two of our cooler targets (KOI-4087 and KOI-4622) 
we gathered near-infrared spectra with the TripleSpec 
instrument on the Apache Point Ob- 
servatory 2.5 m telescope to supplement the material de- 
scribed above and to aid in the determination of the stel- 


Table 4 
Spectroscopic observations (Keck/HIRES). 


BJD Exposure SNR RV 
Star (2,400,0004-) (sec) (pix-!)  (kms-!) 
KOI-0571 56135.11242 900 18 —61.7 
KOI-1422 56192.77304 900 13 —24.8 
KOI-2529 56472.07615 900 13 —12.2 
KOI-3255 56472.06602 729 24 —15.3 
KOI-3284 56450.88292 2700 45 —90.9 
KOI-4005 56472.05630 785 25 —41.0 
KOI-4087 56472.04532 900 15 —20.1 
KOI-4427 56472.03348 900 12 —11.6 
KOI-4622 56472.08746 900 15 —98.6 
KOI-4742 56472.02224 900 20 —76.2 
KOI-4745 56472.09849 900 14 —24.0 


lar properties (Sect. [5). For both observations, taken on 
UT 2014 June 26, we used the 1"1 slit giving a resolv- 
ing power of R ~ 3500. We alternated exposures in an 
‘ABBA’ fashion, integrating for 4 minutes at each po- 
sition until we achieved a SNR of at least 80 per pixel 
throughout the H band (our SNR was slightly higher, 
about 100, for KOI-4622). Our total integrations were 
48 minutes and 64 minutes for KOI-4087 and KOI-4622, 
respectively. We took care to gather spectra of the A0 V 
telluric standard star 51 Dra (using an exposure time of 
10s per image in the same ABBA nod strategy) at in- 
tervals separated by no more than 30 minutes in time 
from the KOI observations, corresponding to changes in 
airmass of 0.1 from the positions of the KOIs. We also 
gathered flat field images using the bright quartz lamp at 
the Apache Point Observatory. We used the spextool 
pipeline to co-add flat fields and 
science frames and to extract the spectra from the full- 
frame images. The xtellcor package was then employed 
to compare the spectra of our standard AO V star 51 Dra 
to that of Vega, in order to identify and remove telluric 
absorption lines in both KOI spectra. 


9. STELLAR PROPERTIES 


We used our Keck/HIRES spectra to estimate the 
effective temperature (Teg), surface gravity (logg), 
and metallicity ([Fe/H]) of our stars using a proce- 
dure referred to as SpecMatch (Petigura et al, in 
prep.). SpecMatch compares an observed target spec- 
trum against a set of approximately 800 library spec- 
tra obtained with the same instrument that span a wide 
range of parameters (Teg = 3500-7500K, logg = 2.0- 
5.0 (cgs), and [Fe/H] values from —1.0 to +0.5 dex). For 
each comparison SpecMatch places the target spectrum 
onto the same wavelength scale as the library spectrum 
and computes x?, the sum of the squares of the pixel-by- 
pixel differences in normalized intensity. We adopted the 
mean Tog, logg, and [Fe/H] of the ten closest matching 
library spectra (weighted by the x?) as the final set of 
stellar parameters for each target spectrum. They are 
reported in Table [5| for about half of our sample. We 
estimate these parameters to be accurate to about 100 K 
in Tog, 0.10dex in logg, and 0.10dex in [Fe/H] based 
on comparisons with benchmark stars having parameters 
measured using detailed LTE modeling. 

For stars cooler than about 4300 K this technique be- 
comes less reliable. In that temperature range the library 
of reference spectra is quite sparse, and it is possible for 
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Table 5 
Spectroscopic parameters. 


Toff log g [Fe/H] 

Star (K) (cms~?) (dex) Source 
KOI-0571 3755+ 90 —0.26 + 0.12 1,2 
KOI-1422 3572+ 80 tt —0.12 + 0.12 1,2 
KOI-2529 4651+100 4.64+0.10 +0.01 + 0.10 3 
KOI-3255 45514100 4.67 + 0.10 0.00 + 0.10 3 
KOI-3284 3748 +112 Z +0.16 + 0.14 2 
KOI-4005 5431+100 4.50+0.10 +0.02 + 0.10 3 
KOI-4087 4134+ 154 Z —0.30 + 0.15 4 
KOI-4427 3813 + 112 —0.07 + 0.14 2 
KOI-4622 4340 + 177 Z —0.57 + 0.18 4 
KOI-4742 4402 +100 4.71 +0.10 0.37 + 0.10 3 
KOI-4745 4723+100 4.621 0.10 0.01 + 0.10 3 


Note. — Sources are: 1. (2013a[B); 2. 
Muirhead et al] (2014); 3. SpecMatch; 4. New estimates from 
this paper. 


the SpecMatch Tig value to be biased high simply be- 
cause the best matches are hotter. Similar biases are 
likely in log g and [Fe/H] Six stars in our sample are in 
this regime: KOI-0571, 1422, 3284, 4087, 4427, and 4622. 
Spectroscopic parameters (except for log g) for the first 
three and for KOI-4427 were reported b 
and/or based on the 
measurement of temperature- and metallicity-sensitive 
features in near infrared H-band or K-band spectra. We 
adopted averages of those determinations here. 

For KOI-4087 and KOI-4622 we derived new param- 
eters based on the NIR spectra described in Sect. [1.4] 
Effective temperatures were inferred from the empiri- 
cal H-band calibrations of [Newton et al] (2015), which 
were established using stars with bolometric luminosi- 
ties and temperatures from long-baseline interferometry 
(Boyajian et al][2013). They are valid between 3100K 
and 4800K, and have a typical scatter of 72 K. Spec- 
tral indices based on the Mg I doublet at 1.50 um 
and the Al I doublet at 1.67 um were measured as de- 
scribed by (2015), and temperature un- 
certainties were obtained from numerical simulations 
following the same work. Metallicities for KOI-4087 
and KOI-4622 were determined using the NIR cali- 
brations of (20134), which have the ad- 


vantage over other calibrations developed in recent 
ears (Rojas-Avala et all[2010, 2012 
2014al that they were established includ- 


ing K dwarfs with effective temperatures greater than 
4000 K. They are thus the most appropriate calibrations 
for these two KOIs, which lie on the K/M-dwarf bound- 
ary. showed that the K-band rela- 
tions display the strongest correlation with [Fe/H]. We 
adopted their calibration in Eq.(16), which has a scatter 
of 0.11dex. The analysis of the spectra, the measure- 
ment of the spectral features (equivalent widths of the 
2.2 014m Na I doublet, the CO (2-0) band-head in the 
K band, and the H20-K2 index from [Rojas-Ayala et al] 
[2013], as well as the uncertainty estimation followed the 


19 Even though SpecMatch has difficulty obtaining accurate pa- 
rameters for cool stars, it is still easily capable of distinguishing 
dwarfs from giants as the pressure-sensitive lines such as the Mg Ib 
triplet are very much narrower in giants than in our dwarf library 
spectra. All of the targets in this paper are verified to be dwarfs. 

20 The [Muirhead et al] (2014) values supersede earlier ones re- 
ported by [Muirhead et al] (2012) that used the same spectra. 


procedures described recently by [Muirhead et al] (2014), 
who applied the same methodology to a similar set of 
spectra of about 100 M-dwarf KOIs. The new tempera- 
tures and metallicities for KOI-4087 and KOI-4622 may 
be found in Table 5] 

The spectroscopic parameters for our targets were used 
to estimate the stellar properties (primarily the mass and 
radius) by appealing to stellar evolution models from the 
Dartmouth series (Dotter et aL]B008). following a Monte 
Carlo procedure similar to that described by Torres ct al] 
(2008). While this is straightforward for the stars with 
complete information (Tog, [Fe/H], and log g), we lack a 
log g estimate for the six cooler stars because the spectro- 
scopic techniques applied to them do not constrain that 
property. Thus, we cannot establish their precise ages, 
or equivalently, their sizes. This uncertainty is relatively 
unimportant, however, because the radii of cool main- 
sequence stars change little with age. Most of those stars 
show periodic brightness variations that may be inter- 
preted as rotational modulation, and in principle knowl- 
edge of the rotation periods (Pot) enables one to infer 
a rough age using gyrochronology relations. The age, in 
turn, may be used in place of logg as a constraint for 
the stellar evolution modeling, although in practice the 
constraint is weak because the radius is not very sensi- 
tive to age, as just stated. Our determination of Pio for 
KOI-3284, 4087, 4427, and 4622 is explained in the fol- 
lowing section; for KOI-0571 and KOI-1422 we adopted 
the periods measured previously by others (see below). 
Our age estimates relied on the gyrochronology relations 
of (2014), which are claimed to 
provide values good to ~1 Gyr for stars with masses 
above 0.55 Mo, degrading rapidly below that. The B-V 
color indices required by these calibrations were taken 
from the U BV survey of the Kepler field by [Everett et al] 
(2013), along with reddening values from the KIC. Gen- 
erous uncertainties of 0.1 mag were assigned to the de- 
reddened color indices to account for possible systematic 
errors. For KOI-1422 and KOI-3284 the colors were ad- 
ditionally corrected for the presence of the close compan- 
ions described in Sect. assuming they are physically 
associated. T'he formal ages inferred for KOI-0571, 1422, 
3284, 4087, and 4622 are 3.89*092, 3.57*325. 4.36* 0.78. 
1.20*025. and 1.82*0 31 Gyr, respectively. As described 
below, KOI-4427 does not display a clear signature of 
rotation, so we have conservatively adopted a broad in- 
terval of possible rotation periods (10—45 days) based on 
the full range of periods observed for M stars of similar 


temperature, as reported by[McQuillan et al] (see 
their Fig. 10). The corresponding age range from the gy- 
rochronology relations of 
is 1-6 Gyr, which we then used as a constraint for the 
isochrone modeling. 

KOI-0571 and KOI-1422 are systems with multiple 
transiting planet candidates (five each). In such cases, if 
the candidates can be assumed to transit the same star 
in near-circular orbits (required for stability), the transit 
light curve modeling can provide a much stronger con- 
straint on the mean stellar density (p,) than for single- 
planet candidates because the density estimates from the 
different candidates within the same system can be aver- 
aged together (see Sect.[8.T]. For these two KOIs we have 
therefore made use of the photometric p, in our Monte 
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Carlo procedure simultaneously with the age constraint 
to strengthen the determination of the stellar character- 
istics. 

The final properties from our stellar evolution mod- 
eling are listed in Table [6] computed from the mode 
of the corresponding posterior distributions. Uncertain- 
ties correspond to the 68.3% (1c) credible intervals from 
the same distributions. Distances were computed from 
the absolute K, magnitudes inferred from the isochrones 
along with the apparent brightness from the T'wo Mi- 
cron All-Sky Survey (2MASS; in the 
same passband, which is the one least affected by in- 
terstellar extinction. We nevertheless made corrections 
for extinction that we inferred from reddening values for 
each KOI adopted from the KIC, as above. For the four 
targets with close companions (KOI-1422, 2529, 3255, 
and 3284) the apparent magnitudes were corrected for 
the light contribution of the neighbors. On the assump- 
tion that those stars are physically bound to the targets, 
their typical orbital semimajor axes would be approxi- 
mately 52, 248, 75, and 64 AU, and their orbital periods 
roughly 550, 4600, 770, and 690 yr, respectively. 


5.1. Stellar rotational periods 


Here we give the details of our determination of the ro- 
tation periods used above for four of the cooler targets in 
our sample that have no spectroscopic estimate of log g: 
KOI-3284, 4087, 4427, and 4622. For this application 
the processing of the raw photometry was somewhat dif- 
ferent than that described earlier, because we wished to 
retain the astrophysical variations present in the light 
curve while at the same time removing the instrumental 
effects. We therefore used the Presearch Data Condi- 
tioning Maximum A-Posteriori (PDC-MAP) data from 
Kepler (Smith et al] 2012), which is designed to meet 
those goals. 

Since the data are unevenly sampled and each quar- 
ter has a unique offset, we elected to use a Lomb-Scargle 
style periodogram[] The light curve model is a simple 
sinusoid and thus is linear with respect to the model 
parameters for any trial rotation period, Pot. Using 
weighted linear least squares we are guaranteed to find 
the global maximum likelihood solution at each trial Pot- 
We scanned in frequency space from twice the cadence 
up to twice the total baseline of observations, taking 10? 
uniform steps in frequency. At each realization, we de- 
fined the “power” as (BICyun— BICtiia)/ BICya, where 
BIC is the Bayesian Information Criterion and “null” 
and "trial" refer to the two models under comparison. 
We also performed a second periodogram analysis with a 
finer grid step around any prominent peaks found in the 
original periodograms. 

The resulting periodograms for KOI-3284, 4087, 4427 
and 4622 are shown in Figure KOI-3284 and KOI- 
4087 show clear uni-modal peaks at 37.53 + 0.73d and 


21 Stars with active regions have a non-uniform surface bright- 
ness distribution, leading to brightness variations as the star ro- 
tates (Buddingl[1977]). These active regions tend to evolve in lo- 
cation and amplitude over timescales of days to years, which can 
cause the periodicities to change as well due to differential rotation 
(Reinhold et al][2013). Despite the complex nature of individual 
spots, the ensemble population tends to imprint the rotation period 
as a dominant peak in the Fourier domain, allowing for an estimate 
of the rotation period using photometry alone (Baer ct al] OTI] 
[Nielsen & Karofil[2012). 
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Figure 4. Periodograms of the Kepler PDC-MAP photometry 
for four of the KOIs for which a spectroscopic surface gravity con- 
straint is lacking. The insets show a close-up of the main peak, with 
an indication of the full width at half maximum that we assign as 
the error in Prot. No clear result is obtained for KOI-4427. 


17.62 + 0.08 d, respectively, which we attribute to their 
rotation periods. Here we have taken the lo confi- 
dence intervals to be the full width of the peak in the 
periodogram at half maximum. For comparison, we 
note that used an autocorrela- 
tion periodogram and found peaks at similar periods 
(37.026 d and 17.542 d). Additionally, KOI-4087 is in the 
sample studied by (3013). who obtained 
Prot = 17.496d, in good agreement with the other re- 
sults. KOI-4427 shows no clear peaks above the noise 
level, in agreement with the findings of 
(2013), so we were unable to assign a rotation period in 
this case. A possible reason for complex structure of the 
periodogram in this case is the presence of other unre- 
solved stars in the aperture. 

KOI-4622 exhibits two peaks and is not included in 
either the or the 
samples. The dominant peak occurs at 21.85 + 
0.22d, with a second peak at ~25% lower power with 
P = 10.72 + 0.06d. We note that this second peak oc- 
curs not only at a lower power but also at nearly exactly 
one half the period of the higher peak. We argue that 
the longer period is that of the stellar rotation, and the 
shorter period is an alias caused by more than one active 
region on the stellar surface. 

Two other cool KOIs have had their rotation periods 
determined previously (33.70 + 0.31d for KOL0571 and 
36.5 + 17.6d for KOI-1422; [McQuillan et al]|2013), and 


we adopt those values here as published. 


6. CANDIDATE VALIDATION 


The Doppler signals (velocity semi-amplitudes) ex- 
pected of the candidates in our sample, if due to planets 
around the targets, range from 0.4 to 1.4ms~! based 
on rough estimates of their masses using the prelim- 
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Table 6 
Stellar properties. 


Age M, Ry, log g 
Star (Gyr) (Mo) (Ro) (cms~?) 
+0.6 +0.024 +0.023 +0.020 
KOI-0571 4.0Ł%6 0.544t9:924  0,523+9:023 . «rag 0.000 
+3.4 +0.033 +0.038 +0.025 
KOl1492 4.2134  0.454t9:033  0.426+%:938 — 4 999-0025 


4.6197 0.015 


+7.7 
KOI-2529 3.0 95 0.028 


+7.5 +0.033 +0.026 +0.015 
KOLe255 2.0752 0.707+%933 Quest? 026 — 4,539*0 015 
T 0.8 +0.041 +0.038 +0.059 
KOL3284 4485 opie  0.520+9:938  4,740+9:059 
+3.6 +0.044 +0.076 +0.035 
KOL4005 7.238  0.884*0044  Q,gggt*O076 — 4514*0035 


KOL4087 1.3799 4.706* 6-16 


+2.6 +0.040 +0.038 +0.067 
KOI-4427 detto pudet) ORI. Avi 
+0.5 +0.049 +0.038 +0.047 
KOr4522 1.995  0572*005  Q5BQ 0055 4715+007 


KOL4742 2.955 
KOL4745 3.2705 


+0.018 
4.673 9:021 
+0.016 
4.614% 9.039 


Px Ly My Mk, Distance 
(gcm-?) (Lo) (mag) (mag) (pc) 
529034 o0 003 — 901024 sarta 172118 
vol oe. (10.12032 6.071929 226175 
305024 0,19949039  g.984025 4.g6#9-11 gists 
318925 0.17349935 7.294938 qato aria 
5.524077  O044' oi, — 9.55 944 — 5.507035 145133 
1897025  0.581+0%153 — mde oo. . dos ois dax eo 
4.767592  0.079*0022.— 8.337033 — 5.02012 — 201*16 
5.797585  0.043*0017 9.344937 — 5.54*032 — 240732 
494*100  Q.0890093 — g.024252 — 5o0*04 2841% 
401 03, OlIT'Dog 7731025 4790m 34213 
206924 02170043 — oasf025 4g2t011 7791% 


Note. 


— Ages, logg, and p, are best-fit values from the Dartmouth models, constrained by either a spectroscopic logg, a 


light-curve-derived px, and/or a gyrochronology age based on the rotation period (see text). 


inary radius measures and the mass-radius relation of 
(2014). Given the faintness of the par- 
ent stars (Kp = 14.3-15.9), the detection of such small 
variations would be very challenging, and therefore the 
planetary nature of these objects cannot presently be 
"confirmed" in the usual way, i.e., by establishing that 
the orbiting objects are of planetary mass. Instead they 
must be “validated” statistically, by showing that the 
likelihood of a true planet (which we refer to hereafter 
as the ‘planet prior’) is orders of magnitude larger than 
that of a false positive. We describe this process below. 


6.1. Overview 


'The types of astrophysical false positives we consider 
in our analysis involve other unseen stars in the photo- 
metric aperture of Kepler that are eclipsed by an orbiting 
object, and have their eclipses attenuated by the light of 
the target such that they mimic shallow planetary tran- 
sits. Examples of these situations include background 
or foreground eclipsing binaries (‘BEB’), background or 
foreground stars transited by a (larger) planet (‘BP’), 
and physically associated stars transited by a smaller 
star or by a planet. Such physically associated stars 
will typically be close enough to the target that they 
are generally unresolved in high-resolution imaging. We 
refer to these hierarchical triple configurations as ‘HTS’ 
or ‘HTP’, depending on whether the object orbiting the 
physical companion is a star or a planet 


The procedure we used for validation is BLENDER 
2011 2013), which 


22 We point out that other statistical studies using Kepler data 
(c.g.,[Morion & Johinsonl[0T T; Morton 2012 [Lissauer et aL[2014) 
have not considered the HTP scenario as a false positive, arguing 
that the exact location of the planet (whether on the intended tar- 
get or a bound companion) is inconsequential although its size may 
be larger if around the companion. Validation without consider- 
ing these kinds of false positives is significantly easier, particularly 
since they tend to dominate the blend frequency, as we show be- 
low. For the present work, however, planet size is critical because 
it affects habitability: a planet that transits an unseen companion 
instead of the target may be too large to be rocky, and is thus less 
interesting for our purposes (see, es. AAA [Rogers[2014]. 


We therefore count HTP configurations as blends. 


has been applied successfully in the past to many 
of the most 


interesting candidates revealed by Ke- 
l 
BLENDER makes full use of the detailed shape 
of the transits to limit the pool of viable blends. It does 
this by simulating large numbers of blend scenarios and 
comparing each of them with the Kepler photometry in 
a x? sense. Fits that give the wrong shape for the tran- 
sit are considered to be ruled out. This enables us to 
place useful constraints on the properties of the objects 
that make up the blend, including their sizes and masses, 
overall color and brightness, the linear distance between 
the background/foreground eclipsing pair and the tar- 
get, and even the eccentricities of the orbits. Those 
constraints are then used to estimate the frequencies of 
blends of different kinds. We note that the simulated 
light curves generated by BLENDER take full account of 
any known extra light in the aperture, such as that com- 
ing from the close companions reported earlier (see Ta- 
ble 2). We also point out that here we only consider 
blends involving main-sequence stars, as BEBs with a 
giant component produce light curves with a shape that 
does not mimic a true planetary transit when observed 
with such high precision as delivered by Kepler. 

The photometric data we use here are the long-cadence 
PDC-MAP time series, as in Sect. [5.1] detrended to 
remove signals at time scales long compared to the 
transit durations. Throughout this analysis we fol- 
lowed the nomenclature established in previous BLENDER 
studies, designating the objects in the eclipsing pair 
as the “secondary” and “tertiary”, and the target it- 
self as the “primary”. We drew stellar properties for 
the primary, secondary, and tertiary (masses, radii, 
and absolute brightness in the Kepler and other pass- 
bands) from model isochrones from the Dartmouth se- 
ries areae Tran the selection of the isochrone 
for the primary was based on the spectroscopic proper- 
ties given earlier. The BLENDER studies cited above may 
be consulted for further technical details. 

In summary, our validations proceed in two stages. 
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First, we use BLENDER proper to derive constraints on 
blend scenarios from their light curve shapes. Then, 
we incorporate other constraints from follow-up obser- 
vations to compute the blend frequencies and the planet 
prior by means of Monte Carlo simulations. We now 
address each of these steps in turn. 


6.2. BLENDER constraints 


To illustrate the constraints derived from the shape of 
the transit, we focus here on one of our candidates, KOI- 
4005.01. BLENDER indicates that background eclipsing 
binaries are only able to produce viable false positives if 
the primary star of the eclipsing pair is restricted to a 
narrow range of main-sequence massed77 between about 
0.7 Mo and 1.4 Mo, and to a limited interval in bright- 
ness (Kp magnitude) relative to the target correspond- 
ing to AKp < 5.5. Figure [5] shows the x? landscape 
for all blends of this kind (BEB scenario) in a repre- 
sentative cross-section of parameter space. Regions out- 
side of the 30 contour correspond to configurations with 
light curves giving a poor fit to the observations, i.e., 
much worse than a true planet fit. These blends are 
therefore excluded. The figure also illustrates some of 
the additional constraints available from our follow-up 
observations for this candidate. For example, analysis 
of our Keck/HIRES spectra of KOI-4005 generally rules 
out companions within 5 magnitudes of the primary if 
they are angularly close enough to fall within the spec- 
trograph slit] This can eliminate much of parameter 
space (see green hatched area in the figur Addi- 
tionally, by comparing the r — K, colors of the simu- 
lated blends with the measured color index of the target 
(r — K, = 1.580 + 0.029;[Brown et al][2011], we find that 
some of the BEB scenarios we have simulated are either 
too blue or too red by more than 3c (blue hatched areas 
in the figure), and are therefore also excluded. How- 
ever, for this particular KOI all BEB blends with the 
wrong color are already excluded by BLENDER for giv- 
ing poor fits (i.e., they are outside of the 3e contour). 
po have pointed out that an addi- 
tional source of blends involves eclipsing binaries with 
eccentric orbits that are oriented so that they show only 
a secondary eclipse, as viewed by Kepler. More gener- 
ally, we note that they could also show only a primary 
eclipse. We find, though, that while the depth and shape 
of the diluted eclipses may indeed match the transit sig- 
nal in some cases, the combined color and brightness of 
these blends are such that they are generally ruled out by 
the spectroscopic constraint and/or the measured color 
of the KOI. 


23 Stars more massive than about 1.4 Mc; will generally have 
evolved to become giants for ages typical of the field (~3 Gyr), and 
such stars yield light curves that have a very different shape than 
a transit, as stated earlier, so they do not constitute viable blends. 

24 For typical seeing conditions at the Keck telescope stars that 
are beyond 0/’43 can still imprint their lines on the target spec- 
trum, and be detected, as can stars at wider separations that hap- 
pen to be aligned along the slit. Nevertheless, to be conservative 
we assume here that these companions are spectroscopically unde- 
tectable. 

25 Note, however, that not all blends in the green hatched area 
are excluded. This is only the case if the background star falls on 
the slit, and its RV is more than 10 km s-! different from that of 
the target. In all other cases we consider the blends to be viable, 
even if they have A Kp < 5. 
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Figure 5. Map of the x? surface (goodness of fit) for KOI- 
4005.01 corresponding to blends involving background eclipsing 
binaries. On the vertical axis we represent the linear distance 
between the BEB and the target (Dpgp — Dtarg), expressed 
for convenience in terms of the difference in distance modulus, 
Aó = 5log( DBEB /Dtarg). Only blends within the solid white con- 
tour (darker colors) provide fits to the Kepler light curve that are 
within acceptable limits (3e, where o is the significance level of 


the x? difference compared to a transiting planet model fit; see 
[Eresarm et al 01D). Other concentric colored areas (lighter colors) 
represent fits that are increasingly worse (4o, 5c, etc.), which we 
consider to be ruled out. The blue cross-hatched areas correspond 
to regions of parameter space where the blends are either too red 
(left) or too blue (right) compared to the measured r — Ks color 
of the target, by more than three times the measurement uncer- 
tainty. The dashed green line labeled A Kp = 5.5 is tangent to the 
white contour from above and corresponds to the faintest viable 
blends. The green line labeled A Kp = 5.0 represents the spectro- 
scopic limit on faint background stars. All simulated blends below 
this line (green hatched region) are brighter and are generally ex- 
cluded if the BEB is angularly close enough to the target to fall 
within the slit of the spectrograph. Thus, very few blends remain 
viable. 


For blends involving a background or foreground star 
transited by a planet (BP scenario) there is a wide range 
of secondary masses that produce good fits to the Kepler 
photometry of KOI-4005, as shown in Figure[6] BLENDER 
indicates that the faintest of these blends are about 5.8 
magnitudes fainter than the target in the Kp passband. 
In this case, however, the spectroscopic and color con- 
straints drastically reduce the pool of viable false posi- 
tives. 

The x? map for blends involving a physically associ- 
ated companion to KOI-4005 transited by a larger planet 
(HTP scenario) is seen in Figure [7] and shows the size 
of the tertiary as a function of the mass of the com- 
panion star (secondary). In this case BLENDER restricts 
the false positives to be in a narrow strip of parame- 
ter space corresponding to secondary masses larger than 
about 0.25 Mo, and planetary sizes between 0.25 and 
about 0.8 Ryup (2.8-9.0 Re). As in the other scenarios, 
color and brightness constraints allow us to reject many 
of these blends. 

Finally, BLENDER indicates that pairs of eclipsing stars 
orbiting the target (HTS scenario) invariably produce 
light curves with the wrong shape for a transit, or fea- 
ture noticeable secondary eclipses that are not observed 
in the photometry of KOI-4005.01, or, if they show only 
a single eclipse due to a high eccentricity and special 
orientation (Santerne et aloa. the overall brightness 
would make the eclipsing binary detectable and/or its 


12 Torres et al. 
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Figure 6. Similar to Figure[5](and with the same color scheme) 
for blends involving background or foreground stars transited by 
a planet (BP scenario). The faintest blends giving acceptable fits 
have A Kp = 5.8 relative to the target (dashed green line). Blends 
below the solid green (AKp = 5.0, hatched green area) are generally 
excluded by the spectroscopic constraint unless the intruding star 
is more than 0"43 from the target, or closer than 0/43 but with 
ARV < 10kms-! (see text). 
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Figure 7. Similar to Figure[5]for the case of physical companions 
to KOI-4005 that are transited by a planet (HTP). Only companion 
stars with masses larger than about 0.85 Mo or smaller than about 
0.42 Mc yield blend colors that are consistent with the measured 
r — Ks index of KOI-4005. However, many of these blends are 
eliminated by the spectroscopic constraint, as they are brighter 
than A Kp = 5.0. 


color inconsistent with the measurement. These types 
of configurations are therefore easily excluded. A simi- 
lar result was found for each of the other candidates in 
our sample, and indeed all previous BLENDER analyses of 
KOIs have also found that HTS blends are always ruled 
out with photometry of the quality delivered by the Ke- 
pler instrument, when combined with the observational 
constraints. Consequently, we do not consider HTS sce- 
narios any further. 

Maps of the x? landscapes for the BEB, BP, and HTP 
scenarios analogous to those described above, along with 
the additional restrictions based on the color and the 
spectroscopic brightness limits, are shown for each of our 
other candidates in Figure set [SHO] 


6.3. Estimating blend frequencies and the planet prior 


The BLENDER constraints just described significantly 
reduce the number of viable BEB, BP, and HTP blends, 
but not all of the ones that remain can be ruled out by 
our follow-up observations. The next step of the analysis 
is therefore to quantify their expected frequencies, which 
we did via Monte Carlo simulations. Specifically, we gen- 
erated large numbers of realistic false positives according 
to each of the blend scenarios, we eliminated those that 
BLENDER tells us would not match the Kepler light curve 
or that would have been detected in our follow-up obser- 
vations, and we counted the survivors. These numerical 
experiments rely on the known distributions of binary 
star properties, the number density of stars near each 
target, and estimates of the rates of occurrence of tran- 
siting planets and of eclipsing binaries from the Kepler 
Mission itself. Those rates of occurrence (as well as any 
dependence they may have on orbital period or other 
properties) are implicit in the lists of KOIs and eclips- 
ing binaries generated by the Kepler team 
that we used in our analysis, 
when normalized by the total number of targets observed 
by the spacecraft. Since we have recently made sig- 
nificant enhancements in these Monte Carlo procedures 
compared to previous applications, we take the opportu- 
nity to describe them in some detail here. 

For the HTP scenario we simulated companion stars to 
each target following the distributions of binary proper- 
ties (mass ratios, orbital periods, eccentricities) proposed 
ty Iaghasa eral] QUA) and sien DUI We 
placed the companions in random orbits around the tar- 
get assuming isotropically distributed inclination angles 
(i) and uniformly distributed longitudes of periastron 
(w), as well as random orbital phases. Other relevant 
properties of these simulated companions (size, absolute 
magnitudes, and colors corresponding to their simulated 
masses) were inferred from the same isochrone used for 
the target, as the two stars are assumed to be coeval in 
this scenario. In particular, we computed the distance 
to the system by imposing the condition that the to- 
tal apparent brightness of the two stars (accounting for 
extinction) must be the same as the Kp magnitude as- 
sociated with the KOI, and this in turn enables us to 
compute the angular separation. We then assigned to 
each of these secondary stars a random transiting planet 
drawn from the actual list of KOIs hosted at the NASA 
Exoplanet Archive (list downloaded on 2014 March 26), 
which ensures that their properties (including any cor- 
relations among them) are as realistic as possible*] We 
accepted only planets with periods similar to the KOI 
under investigation (within a factor of two). The ratio- 
nale for this is that the relevant HTP blend frequency 
is that of configurations that involve planets with peri- 
ods near that of the candidate, as those frequencies de- 
pend strongly on orbital period. We then examined the 
properties of the stellar companions and their planets, 
and rejected scenarios that do not satisfy the BLENDER 
restrictions on companion mass, planetary size, and or- 
bital eccentricity of the planet (related to the transit 


?6 As pointed out by [Lissauer et al] (2014), when attempting 
to validate targets with a single transiting planet candidate it is 
more appropriate to exclude from the list of known KOIs all stars 
with multiple candidates, as this unduly inflates the rates of occur- 
rence of planets (particularly small ones that are more common in 
multiple systems). We have followed this recommendation here. 
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Figure 8. BLENDER x? landscapes for five KOI targets. Each row corresponds to the KOI labeled in the left panel, and shows the three 
blend scenarios. Observational constraints from spectroscopy and color are also shown. See Figures [BH7] for a detailed description. 
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Figure 9. Same as Figure[8]for a second group of five KOIs in the sample. 
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Figure 10. Same as Figure[8]for the remaining KOI in the sample. 


duration), as they yield poor fits to the transit light 
curve. We also rejected configurations with companions 
bright enough that they would have been detected from 
our high-resolution imaging observations, our spectro- 
scopic observations, by our centroid motion analysis, or 
that would yield an overall color for the blend inconsis- 
tent with the measured color of the target 
[2011]. In applying the spectroscopic constraint we ex- 
cluded companions brighter than the limits described in 
Sect. [4.3] only if they are within 0/43 of the target (half 
width of the spectrometer slit), unless their radial veloc- 
ity computed from the simulated orbit around the target 
is within 10 kms-! of that of the primary. In that case, 
blending of the spectral lines of the secondary with those 
of the primary would prevent detection, so we considered 
those blends still viable. Finally, we retained only the 
false positive configurations that would be dynamicall 

stable according to the criterion of 
(1999). Since not every stellar companion will necessar- 
ily have a transiting planet and therefore act as a blend, 
as we assumed above, we adjusted the blend frequencies 
by accounting for the rate of occurrence of planets of each 
size and period as inferred from the KOI list itself. Minor 
corrections were applied to those rates for incompleteness 
and for the incidence of false positives in the KOI list 
by means of simulations as described by [Fressin et all 
(2013). An additional adjustment was made to account 
for the dependence of the transit probability on the size 
of the secondary relative to the average size of the host 
stars in the KOI list. The final blend frequencies were 
obtained by multiplying by the occurrence rate of binary 


stars with planets as a function of binary semimajor axis, 
following (2014). 


'The calculation of the frequency of background or fore- 
ground stars transited by a planet (BP) proceeded in 
a similar fashion, and scales with the number density 
of stars in the vicinity of each target, for which we re- 
lied on predictions from the Galactic structure model of 


(2003 (Besancon model)P7] We gen- 


27 We note that the current online version of the Besangon simu- 
lator (available at http: //model.obs-besancon.fr ) relies on stellar 
evolution models by |Haywood| (1994) to generate some of the stel- 
lar properties returned, including absolute magnitudes, colors, tem- 
peratures, and radii. Many of the simulated main-sequence stars of 
interest for our work are significantly cooler than the Sun (spectral 
type K and M), and the [Haywood models are not specifi- 
cally designed for such stars as they adopt gray boundary condi- 
tions and a rather simple equation of state, both of which result in 
increasingly biased predictions for cool stars (e.g., 
[2001]. Because of this, we have preferred to recompute the relevant 
properties of all simulated stars adopting from the output of the 
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erated a list of simulated stars in a 5 square-degree area 
around each of our targets, including their kinematic 
properties (radial velocity), and drew stars randomly 
from this list assigning them a random angular separa- 
tion from the target within the corresponding 3c exclu- 
sion region derived from our centroid motion analysis 
(since stars outside this area would have been detected; 
see Sect. [4.2). We assigned a random planet to each of 
these stars that we drew from the KOI list, keeping only 
those with orbital periods within a factor of two of that of 
the candidate, as done previously. False positive scenar- 
ios that do not meet the constraints from BLENDER were 
rejected, along with those that would have been flagged 
by our high-resolution imaging observations or color in- 
formation. Further blends were excluded by the spectro- 
scopic limits on the brightness of any unseen companions, 
as before, except for background/foreground stars whose 
spectral lines would be blended with those of the pri- 
mary (i.e., those with simulated radial velocities within 
lOkms-! of the measured heliocentric velocity of the 
target), which we regarded as viable. Surviving blends 
were weighted by the corresponding rate of occurrence of 
transiting planets of that size and period, corrected as in 
the HTP case for incompleteness and false positives and 
adjusted for the dependence of the transit probability on 
the secondary radius. The final blend frequencies were 
then normalized by the ratio of the areas between the 
centroid exclusion region and 5 square degrees. 

'To estimate the frequency of BEBs acting as blends 
we again drew stars randomly from a list generated 
with the Besançon Galactic structure model, and paired 
them with an orbiting stellar companion (tertiary) gen- 
erated from the distributions of binary properties by 
and [Tokovinir] (2014). We as- 
signed orbital periods to these secondary-tertiary pairs 
drawn randomly from the catalog of Kepler eclipsing 
binaries of (2011), retaining only those 
within a factor of two of the period of the candidate. 
Each blend was also assigned a random angular separa- 
tion from the target within the 3c exclusion limit from 
our centroid motion analysis. After rejecting scenarios 
that do not meet the constraints from BLENDER for ac- 
ceptable fits to the transit light curve, we applied the 


Besangon simulator only the basic characteristics that derive from 
the initial mass function and the Galactic structure model, which 
are the stellar mass, distance, age, and metallicity. We then used 
model isochrones from the Dartmouth series (Dotter ct al [2008) to 
generate all other stellar quantities for our purposes, as the phys- 


ical ingredients of these isochrones give more realistic predictions 
for the lower main-sequence. 
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observational constraints (imaging, spectroscopy, color 
information) as in the BP scenario, and tallied the re- 
maining viable blends. For computing the final BEB fre- 
quency, each surviving blend was assigned a weight given 
by the ratio of the eclipse probability of the particular 
blend (a function of the secondary and p En 
to the average eclipse probability for the |Slawson et al. 
sample as a whole. 

'The probability that the candidate is a true planet as 
opposed to a false positive (‘planet prior’, PL) was es- 
timated by simply counting the number of known KOIs 
with periods within a factor of two of each candidate, 
radii within 30 of the measured planetary radius, and 
with similar transit durations. Incompleteness and false 
positive corrections were applied in the same way as 
above. 


6.4. Results 


Blend frequencies for the HTP, BP, and BEB scenar- 
ios are presented for each of our candidates in Table 
along with the planet priors, PL. For the present work 
we have set a threshold for validation equivalent to a 3a 
confidence level, consistent with previous applications of 
BLENDER. The statistical significance may be expressed 
as PL/(PL + HTP + BP + BEB). For a 30 valida- 
tion we therefore require a planet prior that is at least 
1/(1/99.73% — 1) zz 370 times larger than the total blend 
frequency. We list these ‘odds ratios’ in Table [7| with 
the corresponding significance levels. For eleven of the 
twelve candidates we have achieved very robust valida- 
tions at the 3c level or higher, strongly supporting their 
true planetary nature. On this basis, the eight that have 
not been previously validated by others are given the new 
planet designations Kepler-436 b to Kepler-443 b (see Ta- 
ble[I). KOI-4427.01 has a lower odds ratio of only ~100, 
corresponding to a significance level slightly over 9996. 
To avoid confusion in the rest of the paper we will con- 
tinue to refer to the newly validated planets by their KOI 
numbers. 

Strictly speaking, the above results compare the likeli- 
hood that the transit signals originate on the target itself 
with the likelihood that they come from eclipses asso- 
ciated with some other unseen star in the photometric 
aperture. However, four of our target stars (KOI-1422, 
2529, 3255, and 3284) have known close companions from 
our high-resolution imaging observations, none of which 
can be ruled out as the source of the transit signal by our 
centroid motion analysis in Sect. Thus, the possibil- 
ity remains that the planets causing the signals orbit the 
companions rather than the targets. In that case, the 
inferred radius of the planets would be somewhat larger 
than the values derived in the next section, the exact 
amount depending on the brightness difference (dilution 
effect) and the physical size of the companion stars. 

An estimate of the likelihood that the planets tran- 
sit the companions was derived as follows. In each of 
the four cases we made the assumption that the com- 
panions are physically associated with the target. This 
is reasonable, as the probability of having an unrelated 
background star as close to the primary stars as observed 
is very small compared to that of a physical companion 
(see also[Horch et_al./2014), based on the number density 


of stars around each KOI from the Besancon model and 
the known rate of occurrence of binaries (Raghavan et al. 


[2010 [Wang et alJ[2014]. Furthermore, no signs of any of 
these companions are seen in our Keck/HIRES spectra 
even though in most cases they are bright enough to have 
been detected. This implies that their radial velocities 
must be similar to those of the primary stars so that the 
spectral lines are blended, preventing detection. This 
again argues for physical association|^] In each case we 
then used the same isochrone as for the primaries to infer 
the companion properties, and we computed the 'planet 
prior’ in the same way as for the targets, but this time 
assuming the planet transits the companion (i.e., increas- 
ing the dilution factor by the appropriate amount). We 
represent this planet prior for the companions as PLeomp, 
and we report these values in Table[7] Also given in the 
table are the probabilities that the planets transit the 
target rather than the companion, which we computed 
as P[targ| = PL/(PL + PLeomp). For KOI-2529.02 the 
results allow us to state with high confidence that the 
planet orbits the target. For KOI-1422.04, KOI-1422.05, 
and KOI-3284.01 the calculations indicate a slight or 
modest preference for the planet being around the tar- 
get rather than the companion, whereas for KOI-3255.01 
the result is inconclusive. We return to this issue later 


in Sect. [8.3] 


7. LIGHT CURVE FITS AND PLANETARY PARAMETERS 


Having validated eleven of our candidates as bona-fide 
planets with very high confidence, and KOI-4427.01 to 
a lower degree of confidence, we subjected their light 
curves to a detailed modeling to infer the properties 
of the transiting objects. For this we used the de- 
trended SAP time series from Kepler described earlier 
in Sect. The transits were modeled using the stan- 
dard TAT algorithm employing the 
quadratic limb-darkening law. This simple model as- 
sumes a spherical, opaque planet transiting a spherically 
symmetric luminous star on a circular Keplerian orbit. 
We re-sampled the long-cadence data onto short-cadence 
sampling following the method described by 
(2010). to avoid smearing effects. Our model has eight 
free parameters in total. These are the orbital period, P, 
the time of transit center, 7, the planet-to-star radius ra- 
tio, R,/R,, the mean stellar density, p, phot; the impact 
parameter, b, the logarithm of fractional light contamina- 
tion, log 8, and the quadratic limb-darkening coefficients 
qı and q2. All of these parameters have uniform priors in 
our fits, except Px,phot, for which we employed a Jeffrey's 
prior between 10^? and 10? gcm ?, and log 8, for which 
we used an informative Gaussian prior where appropri- 
ate. 

In our model, 8 is the flux of any contaminating sources 
in the aperture divided by the flux of the target star. The 
purpose of this term is to account for dilution of the tran- 
sit light curve due to close companions identified by the 
high spatial resolution imaging observations. We com- 
puted these ratios here based on the magnitude differ- 
ences of the companions listed in Table[2] transformed to 
the Kp band, and assuming that all their flux is included 
in the aperture. For the wider companions (4—5" separa- 
tions) that are also in the aperture the transformed mag- 


?8 Tn an independent study of KOI-1422, (2014) also 
reached the conclusion that the companion is physically bound, 
based on numerical experiments with model isochrones. 
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Table 7 
Blend frequencies, planet priors, odds ratios, and significance level of the validation for our targets. 


Candidate HTP BP BEB PL Odds ratio Significance PLeomp "P[targ] 
KOIL0571.05 1.15 x 10-6 1.32 x 1077 6.60 x 10-8 2.15 x 10-3 1595 99.9496 ee ee 
KOI-1422.04 7.75 x 1077 4.74 x 10-8 1.28 x 10-8 3.19 x 1074 382 99.7496 5.94 x 10-4 78.4% 
KOI-1422.05 4.49 x 1077 4.84 x 10-8 1.29 x 10-8 4.63 x 1074 907 99.89% 7.05 x10-5 86.5% 
KOI-2529.02 3.24 x 1076 2.17 x 10-? 5.55 x 1079 1.31 x 1073 403 99.75% 5.61 x 10-96 99.6% 
KOI-3255.01 4.29 x 1077 4.79 x 10-11 en 4.91 x 107-4 1144 99.91% 4.39 x 1074 52.8% 
KOI-3284.01 6.71 x 1077 1.10 x 107° 9.27 x 10-19 1.57 x 1073 2333 99.96% 8.02 x 10-4 66.2% 
KOIL4005.01 3.71 x 1077 e en 6.61 x 10-4 1782 99.9496 ee e 
KOI-4087.01 1.34 x 1077 s ees 1.32 x 10-4 985 99.90% 

KOI-4427.01 2.23 x 10-8 1.07 x 10-9 1.39 x 10710 2.62 x 10-4 117 99.16% 
KOI-4622.01 9.98 x 10-8 3.82 x 10710 6.68 x 1071? 2.08 x 10-3 20761 99.99% 
KOI-4742.01 1.37 x 10-6 4.89 x 10-19 5.45 x 1071? 1.26 x 1073 919 99.89% 
KOLA4745.001 1.53 x 10-9 1.41 x 10-9 9.19 x 10-7? 2.85 x 1073 1861 99.95% 
Note. — For several of the candidates the blend frequencies for the BP and/or BEB scenarios are negligible and are not listed. 


For the four targets with close companions PLeomp is the planet prior assuming the planet transits the close companion, and P[targ] 
the probability that the planet transits the target rather than the companion. 


nitudes were taken from the CFOP Web site, and for the 
closer ones we carried out the conversion using the Dart- 
mouth isochrone for each target, on the assumption that 
the companions are physically associated. In the case 
of KOI-1422 we adopted the magnitude difference in the 
Kp band as reported M EET (2014). 'The 8 factors 
range from 0.047 + 0.009 for KOI-2529 to 0.83 + 0.18 for 
KOI-3255, which has a very bright companion. Addition- 
ally, other nearby Kepler sources may contaminate the 
light curve (e.g., bleeding into the photometric aperture) 
and were corrected for using the reported contamination 
factors from the Kepler pipeline, treating them as fixed 
parameters. 

It is worth noting that we did not directly fit for the 
standard quadratic limb darkening coefficients uj and 
uz, but instead we used the transformed parameters qi 
and q2 as advocated by [Kipping] (2013b), in order to im- 
pose efficient, uninformative, and physical priors for the 
limb darkening profile. We note also that the out-of- 
transit baseline flux for each transit epoch was fitted as 
well. However, in this case we used a linear minimiza- 
tion to derive the baseline flux, similar to that described 
by[Kundurthy et al] (2011). 'This treats the baseline flux 
simply as a nuisance parameter that is not marginalized 
against, but rather minimized at each Monte Carlo re- 
alization, reducing the number of free parameters yet at 
the same time allowing us to account for any residual 
offsets in the detrended light curves. 

'To regress our 8-parameter model to the observations, 
we employed the multi-modal nested sampling algorithm 
MULTINET described by add 
(2009). We used 4000 live points with con- 
stant efficiency mode turned off, and set an enlargement 
factor of 0.1. The maximum a posteriori model param- 
eters and their associated 68.396 credible intervals are 
collected in Tables [8] and Q] Also found there are other 
derived parameters including the scaled semimajor axis 
(a/R,), the orbital inclination angle i, the transit du- 
rations T14 and Tos (first-to-fourth and second-to-third 
contacts), and the planetary radius and semimajor axis, 
Rp and ap. Additional derived quantities of interest are 
discussed below. The folded transit light curves are dis- 
played graphically in Figure [LJ] along with the models. 

The planetary radii for our targets are all under 2.4 Rg 


with the exception of KOI-2529.02, and the semimajor 
axes range from about 0.15 to 0.64 AU. We find that in 
all cases our planetary radii are larger than those cur- 
rently held on the NASA Exoplanet Archive. This may 
be partly due to the fact that our stellar radii reported 
earlier also tend to be larger, and that in four cases we 
have included dilution effects from close companions. 
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Table 8 


Final parameter estimates for the objects studied in this work. 
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+0.0013 
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--0.0036 
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+0.171 
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+0.022 
0.734 0.020 


KOI-1422.04 


+0.00050 
63.33600 10.00050 
+0.0026 
55692.3051 100026 
+0.0022 
0.03627 0.0018 


+0.29 
1.087 9/15 


KOI-1422.05 


Fitted parameters 


+0.0029 
0.0297 20.0037 
OTE Oa 


+0.26 +0.30 
0.09 0.09 0.182018 
+0.16 +0.07 
0-03 903 0.93 9:18 
+0.22 +0.03 
0.032003 0.97915 


(—0.532 + 0.089) 


Other transit parameters 


137434 Tp 
0.05005 9.890. 
0.3540 Ls2*032 
—0.017 5:30 0.70" 0:37 
3300049 3.62+9:26 
$906 048 3.341931 


—0.041 


Physical parameters 
+0.16 
1.48 10.25 
+0.035 
0.149" 5024 
+0.069 
0.915 0.064 


intr reo Seats id auch de 0.04 007 9127010 0.1079 +5 
Seg [oe | cux stets 0.3079 38 0.3970 78 1.4576-39 
PHZ] 2222225 98.4 99.7 82.0 
[rocky] [A] .........sssse 68.4 30.6 50.7 
Note. — All quantities correspond to the mode of the posterior distributions. 


* Mean stellar density derived from the light curve fit. 


b Gaussian prior. 


* Mean stellar density derived from our stellar evolution modeling. 
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Table 9 


Final parameter estimates for the objects studied in this work. 


KOI-4005.01 


+0.0016 
178.1396 7 9 6018 
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Plrocky] [%] ............0. 6.5 29.8 27.3 
Note. — All quantities correspond to the mode of the posterior distributions. 


* Mean stellar density derived from the light curve fit. 


b Gaussian prior. 


* Mean stellar density derived from our stellar evolution modeling. 
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8. CONSTRAINTS FROM ASTERODENSITY PROFILING 


A comparison between the mean stellar density in- 
ferred from our light-curve fits above (photometric den- 
sity Px,phot) and some independent measure of the den- 
sity, such as that derived from model isochrones, can pro- 
vide very useful information about a planetary system. 
This technique, referred to as Asterodensity Profilin 
(AP), was discussed in detail by[Kipping et al] (2012) 4 
For a “vanilla” exoplanet (circular Keplerian orbit, no 
blends, etc.) these two density metrics should agree. De- 
viations between them occur when one of the following 
physical effects is not accounted for in the light curve 
analysis: i) orbital eccentricity (‘photo-eccentric effect’); 
ii) contaminating light (‘photo-blend effect’); iii) unoc- 
culted star spots (‘photo-spot effect’); iv) transit tim- 
ing variations (‘photo-timing effect’); v) transit duration 
variations (‘photo-duration effect’); vi) substantial mass 
of the transiting object (‘photo-mass effect’); or vii) the 
object orbits a different star. These effects have been 
described in detail by (2014al), who provided 
approximate expressions for each. 

In what follows we apply AP to 1) test the hypothesis 
that all planets in the multiple transiting planet systems 
KOI-0571 and KOI-1422 orbit the same star (Multi-body 
Asterodensity Profiling, or MAP; [2012); 
2) derive minimum eccentricities for the twelve KOIs in 
our sample; and 3) evaluate the possibility that the val- 
idated planet orbits the close companion instead of the 
target, in the four cases where high-resolution imaging 
has identified such sources. 


8.1. Using AP to derive average stellar densities for the 
multi-planet systems KOI-0571 and KOI-1422 


Of the twelve candidates studied in this work, three 
reside in systems with five transiting planets each: one 
orbits KOI-0571, and two are associated with KOI-1422. 
These rare systems provide extra information that we 
will exploit here. For each planetary candidate in each 
system we performed an independent detrending and 
subsequent light curve fit using the methods described 
earlier. 

If all five planets in each of these two systems orbit the 
same star, then we expect all of the planets to maintain 
low orbital eccentricities in order for these compact sys- 
tems to be dynamically stable over Gyr timescales (see 
also[Rowe et al][2014). Therefore, under the assumption 
that the planets orbit the same star, one would expect 
negligible photo-eccentric effects, implying that the mean 
stellar densities derived from the light curves should be 
consistent. On the same assumption that all five planets 
orbit the same star, two other AP effects listed above 
can also be neglected: the photo-blend effect, and effects 
from the planets orbiting different stars. The latter sce- 
nario is obviously not relevant under the posed assump- 
tion, and the former can be negated since, to first order, 
contaminating light that is unaccounted for disturbs the 
light-curve-derived stellar density for all five planets to 
the same degree. 

A comparison of the photometric stellar densities for 
the five planets around KOI-0571 and KOI-1422 reveals 


29 A similar concept was described by [Tingley et al] (2011]), al- 


though the current implementation is considerably more advanced. 


excellent agreement in both cases, as shown in Figure[12] 
Furthermore, similar diagrams (not shown) comparing 
the limb-darkening parameters u; and us also show very 
good agreement among all planets in both systems. We 
conclude from this that each system is consistent with 
having five planets orbiting the same star, given the avail- 
able data. 

Averaging over the four other planets in KOI-0571, 
and the three others in KOI-1422, we may derive a 
mean photometric density for the host stars that we de- 


note px, MAP. We obtained p4, MAP = S4 a gom * 


and 7.2129 gcm™?, respectively, where the extra dilution 
from the close companion to KOI-1422 has been properly 
accounted for, as before. The lo and 2ø confidence inter- 
vals of these values are marked in Figure[12] with dashed 
lines. Both of these mean densities are consistent with 
the corresponding photometric densities derived from 
the light curves of KOI-0571.05, KOI-1422.04, and KOI- 


1422.05 individually, which are o, phot = 5.9* 22 g cm3, 


11.4* T5 gcm’, and 5.321 gem, respectively. 

'The average stellar densities for these multi-planet sys- 
tems (computed from the remaining planets in each case, 
excluding the validated planets) were used as luminosity 
indicators to strengthen the determination of the host 
star parameters via our stellar evolution modeling in 
Sect. [5] (Sozzetti et al][2007). This is particularly use- 
ful for these three KOIs as neither of the host stars has 
a spectroscopic determination of log g available to other- 
wise constrain the luminosity (see Table [5). 


8.2. Using AP to measure minimum eccentricities 


A measure of the stellar density that is independent 
of that derived from the individual light curves of the 
validated planets is available for each of the twelve 
KOIs studied in this work (see Sect. 5). For KOI-0571 
and KOI-1422, this comes from the MAP-based density 
(px MaP) refined by our isochrone analysis by making 
use of additional constraints on age from gyrochronology, 
along with the spectroscopic temperature and metallic- 
ity estimates. For KOI-3284, 4087, 4427, and 4622 our 
stellar evolution modeling used only the age along with 
Teg and [Fe/H], and for the remaining targets that are 
earlier in spectral type we used logg as the luminosity 
indicator, which for these stars we were able to determine 
spectroscopically using SpecMatch. 

'These independent stellar densities, denoted here sim- 
ply by p,, may be compared with the light-curve-derived 
stellar density of each validated planet, px, phot- In these 
cases we expect the photo-eccentric effect to be the dom- 
inant AP effect since contamination has been accounted 
for in the fits (photo-blend effect), the stars show rel- 
atively low activity (photo-spot effect), there are no 
known timing effects for these objects (photo-timing and 
photo-duration effects), and their small sizes suggest they 
should be of low mass (photo-mass effect). We therefore 
attribute any differences between px, phot and p, as being 
due to orbital eccentricity. The minimum eccentricity 
for each KOI is directly given by this comparison via the 
expression presented Mrene (2014al) , 

1 — (Px,phot /Px)?/3| 
1 + (Px,phot/Px)?/? 


Cmin = 
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Figure 11. Light curve fits and Kepler observations (small dots) for our twelve KOIs. Pink circles represent 10-point binned data. 


[8] and [9] 


reveal that most of the KOIs in the sample are consistent 
with orbiting their parent stars on low-eccentricity orbits. 


'The results of this calculation, shown in Tables 


8.3. Using AP to test blend scenarios 


For the four KOIs found to have close companions from 
our high-resolution imaging (KOI-1422, 2529, 3255, and 
3284) the same techniques from the preceding section 
enable us to revisit the possibility that the planets orbit 
the companions rather than the primary stars. As done 
earlier, we made the reasonable assumption that these 
close neighbors are bound to the targets. This allows us 
to reuse the target isochrones to infer approximate stel- 
lar properties for the companions, including their mean 
densities, on the basis of their brightness differences com- 
pared to the target. We re-fitted the transits adopting re- 
vised blend factors (3) appropriate for the different host 
stars to derive photometric densities for the companions, 
Px,phot. We then used these densities to estimate a min- 
imum eccentricity for the scenarios in which the planets 


orbit the companions rather than the brighter stars. 

We report the results of this calculation in Table [10] 
where we list also the planetary radii that we infer if 
the planets transited the companions, as well as their ef- 
fective insolation level (see next section). Although the 
eccentricity distribution of small planets around M dwarf 
stars is presently unknown (Kipping 201). in general 
we expect circular orbits to be more likely than eccen- 
tric orbits, Qe on evidence from other sub-populations 
of planets ( see, e.g., [Kipping||2013a). We may thus use 
€min aS a rough criterion to decide which location for 
the planet is favored. Additionally, the sign of the quan- 
tity log(p. pnot/p«) conveys useful information. From ge- 
ometric arguments, a planet in an eccentric orbit has 
a higher p un of transiting near periapsis than 
apoapsis ( n For the more probable case 
of O<w » ie “near” periapsis transits), one expects 
log(p.., phot / Px) d 0. In contrast, for the less probable 
scenario of * < w < < "near" apoapsis transits), one 
expects log(p, pnot/p«) < P (see [Kipping|[20144]) ps 
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Table 10 
AP-derived minimum eccentricities and other properties for 
KOIs with close companions. 


Parameter Primary Companion 
KOL-1422.04 
Kp [mag].............. 16.15 + 0.03 17.72 £0.15 
a MM (—0.628 + 0.061) (0.628 + 0.061) 
log(p. [£cm-?]) ....... UN I cS 1342400 
log(Px,phot [gem™*)) .. 1.081973 1.084575 
log(Px,phot/Px) -+--+ 0.164077 —0.271 54 
Lm O12 70 ae 0:20 10:20 
Re [Rep] RN rue 2.01405 
Bap [8] festo eges 028") i. UP T 96 
KOI-1422.05 
Kp [mag].............. 16.15 + 0.03 17.72 £0.15 
T1 TR (—0.628 + 0.061) (0.628 + 0.061) 
log(p. [E£cm-?]) ....... 0.915") 054 1342" E 
log(p. pho: [Ecm ?]) .. 0.77019 0.75038 
log(Px,phot/Px) -+--+ —0.14* 03 —0.56*527 
inne sa e De See RISE RS 0.104915 Ug S 
Ry Db receta 142ta 1.68+0:27 
Sag [S5] exire 125" 5 oe 101759 
KOI-2529.02 
Kp [mag].............. 15.86 + 0.10 22.00 + 0.20 
MOE TNNT (—1.330 + 0.084) — (2.474 + 0.089) 
log(p« [gcm-?]) ....... 0 4820038 T 866" 0-040 
log(Ps,phot [g cm ?]) .. 0.72* 055 OBE NTs 
log(Px, pnot/p«) ----..-- 0.24+0:39 -1.1*012 
E Diac EUM 019 0.65+0:06 
Ry E EE PREM ares 
Aust (Be eese nns 1.697 5 75 0.304906 
KOI-3255.01 
Kp [mag].............. 15.01 + 0.10 15.21 + 0.20 
ISP codernecees ata (—0.080 + 0.089) (0.080 + 0.089) 
log(ps [gcm-?]) ....... DIETE Qt ue 
log(p. phot [pcm ?]) .. 0.50* 0:15 0.51* 918 
logos stat DA) eo eres 0.014022 0.0040: 1 
On dee pera ewes ute 0.02+0:03 0:007 506 
Rp [Ra erecto rre ots gout 
Seg [S osea 219. 151704) 
KOL 3284.01 
Kp [mag]... cued 14.58 + 0.10 17.01 + 0.20 
T MENS (—0.577 + 0.070) — (1.036 + 0.088) 
log(ps [gem73])....... QT1T Dore nore] 
log(Px,phot [g m-?]) .. 0.78* 528 0.77 F638 
log(Px, phot/P») -<---> —0.03+9:31 —0.541 0-52 
CH 0.08775 0.407922 
Rp [Ra] tete Lie pgs 
Boe [S0] ones 144p 97 0.864039 
Note. — The log 8 values indicated in parentheses were used as 


Gaussian priors. For KOI-2529.02 and KOI-3284.01 the 8 factors 
for the companions are not the reciprocal of the values for the 
primaries because of the presence of additional companions in the 
aperture (see Table[2). Those additional stars are ruled out as the 
source of the transit signals by our centroid analysis of Sect. 
but still cause dilution. The symbol p, represents the mean stellar 
density based on our isochrone analysis (and other constraints), 
while p, phot is the photometric density from the individual light 
curve fit for the planet. 
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Figure 12. Mean stellar densities from separate light curve fits 
to each of the planets in the five-planet systems KOI-0571 and 
1422. The dashed lines and gray areas represent the lo and 26 
confidence regions of the average density computed from the other 
companions in each system (open circles) excluding the ones we 
have validated here (filled circles). In both cases the agreement 
between this mean density and the separately determined one of 
the validated planets is excellent, supporting the notion that in 
each system all five planets orbit the same star. 


fore, in general the photo-eccentric effect is more likely 
to (but will not necessarily) produce a positive sign for 
this logarithmic ratio. Additionally, a diluted transit 
in a circular orbit will cause log(p,.55ot/p.) to be nega- 
tive via the photo-blend effect eee SEE) but can 
never yield a positive sign. For these reasons, finding 
log(Px,phot/P*) > 0 is more consistent with a bona-fide 
planet than a negative sign. 

As seen in Table[I0] despite the large uncertainties the 
minimum eccentricities are considerably higher if KOI- 
1422.05, KOI-2529.02 and KOI-3284.01 orbit the sec- 
ondary star, which seems a priori less likely, although 
quantifying an accurate odds ratio is not possible without 
knowing the eccentricity distribution of terrestrial-sized 
planets around M dwarfs. For the first two cases this 
assessment agrees with our conclusion from the BLENDER 
study (Table[7), which also favored the planet transiting 
the primary. For KOI-2529.02 and KOI-3284.01 the sign 
of the log-density ratio log(p.pnot/p.) is also negative 
if the planet is around the companion, requiring either 
geometrically disfavored near-apoapsis transiting plan- 
ets, or simply the photo-blend effect. For KOI-1422.04 
the minimum eccentricities between the two scenarios are 
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comparable, but the sign of the log-density ratio favors 
the primary star as the host. Finally, the test for KOI- 
3255.01 is ambiguous according to both metrics, in agree- 
ment with the indications from BLENDER that were also 
inconclusive for this planet. 


9. HABITABILITY 


While a broadband transit light curve alone does not 
reveal the “habitability” of a given planet, it does allow 
us to determine how close to the “habitable zone” a par- 
ticular planet is. Therefore, an important distinction to 
make is that we do not claim any of the planets in our 
sample are truly habitable, but merely that the insola- 
tion they receive from the host star is suitable for the 
presence of liquid water on the surface across the range 
of plausible atmospheric conditions. 

have presented calculations of 
the habitable zone for a wide variety of atmospheric con- 
ditions, and considered also empirical boundaries for the 
hot inner edge modeled by a recent- Venus-like world, and 
a cold outer edge modeled by an early-Mars-like world. 
More PME, ecaer al (2013) took a broader view of 
the conditions that might provide habitable surface con- 
ditions closer to the host star, and postulated an inner 
boundary for the HZ that can be as close as 0.38 AU for 
a solar-type star if the relative humidity of the plane- 
tary atmosphere is low and the albedo high. Given the 
complexity of the problem and our still limited state of 
knowledge, in this work we have chosen to adopt the 
broadest possible limits that can plausibly lead to suit- 
able conditions for life, in order to avoid prematurely 
dismissing as uninteresting planets that may yet be hab- 
itable. Specifically, we adopted the empirical cold outer 
boundary for the HZ of modeled 
by an early-Mars-like world, and the inner edge for dry 
desert worlds from|Zsom et al] (2013). 'These boundaries 
depend upon two key input parameters: Seg, the effec- 
tive insolation received by the planet, and the spectral 
type of the star, for which the usual proxy is Teg, the 
star’s effective temperature. We define Seg in terms of 
the Earth's insolation, Sg, as 


Seg m L,/Lo 
Ss  (a,/AU?2V1—e& ' 


where L, is the stellar luminosity and a, the planetary 
semimajor axis in astronomical units. For each planet 
in our sample we computed a joint posterior for Seg and 
Terr based on the stellar properties inferred in Sect. [5] as- 
suming the planet/s orbit is circular. At each realization 
in this posterior we evaluated whether the sample falls 
inside or outside the adopted HZ boundaries. Count- 
ing the total number of cases inside allows us to assign 
a statistical confidence that the planet in question is in 
the HZ, denoted P[HZ]. We provide those results in Ta- 
bles [8] and [9] This statistical approach to the habitable 
zone was first described by for the 
planet Kepler-22 b. 

Additionally, we consider the radius of the planet 
to be an important parameter in assessing its poten- 
tial for habitability, as this can determine whether 
the planet is likely to have a solid surface that 
can support liquid water. In the past it was of- 
ten conventionally assumed that there likely existed 
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Figure 13. Joint posterior distribution of the planetary radius 
(y axis) and effective insolation (x axis) for KOI-0571.05 (Kepler- 
186f). Labeled vertical lines mark the various boundaries of 
the habitable zone as defined by and 
(2013). Here we adopt the optimistic ‘dry desert’ and 
'early Mars’ inner and outer edges, respectively. These lines are 
drawn assuming an effective stellar temperature equal to the mode 
of the associated posteriors. Contours enclose 68.396, 95.496 and 
99.7% of the Monte Carlo samples (corresponding to the lo, 26, 
and 3c levels). 


some boundary in radius between rocky planets and 
mini-Neptunes (e.g. , perhaps around 
1.75 Rẹ (Lopez & Fortney . However, recent 
counter-examples to this paradigm make such a sim- 
plistic division questionable. For example, KOI-314c 
is a mini-Neptune despite being 1.6 Ra 
[2014al), and Kepler-10 c has been claimed to be rocky de- 
spite being 2.35 Ra (Dumusque et_al./2014; for a differing 
view see and|Wolfgang & Lopez[2014). Nev- 
ertheless, it is reasonable to assume that smaller planets 
are more likely to be rocky, and thus we present joint pos- 
terior distributions in radius and insolation. An example 
of this diagram for KOI-0571.05 is given in Figure[I3] in 
which we indicate various boundaries for the HZ as de- 
fined by[Kopparapu et al] and[Zsom et al] (2013). 
Similar graphs for the remaining targets in our sample 
are shown in Figure [14] 

From the standpoint of habitability it is of interest 
also to quantify the likelihood that the planets are rocky, 
given their measured sizes and corresponding uncertain- 
ties. To do this we have considered an empirical den- 
sity model giving a statistical estimate of the mass of 
a planet as a function of its radius, constructed based 
on the available sample of transiting planets with mea- 
sured densities (i.e., with masses measured either via the 
Doppler technique, or from TTVs). This model is ex- 
pected to be biased towards larger and more massive 
planets that are easier to detect, but it is possible to re- 
move the bias by convolving the density model with es- 
timates of the rate of occurrence of planets as a function 
of their size (Fressin et a1]|2013), which we have done in 
a Monte Carlo fashion (Fressin et al]2015). In this way 
we may infer the true distribution of planets in the mass- 
radius diagram, which can then be used to establish the 
fraction of rocky or gaseous planets at a given radius or 
mass. For the purposes of this work we have considered 
any planet denser than a composition of 100% MgSiO3 
perovskite to be rocky (Seager et al][2007]). Using the 
above model and our posterior distributions for the ra- 
dius of the planets in our sample, we have computed the 
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Figure 14. Similar to Fig. [[3]for the other planets in this study. 


probability that they are rocky, P[rocky], by assigning a 
mass to each link in the chain and counting the number 
of samples denser than the threshold, analogous to what 
was done previously for the HZ. The P [rocky] values are 
2 for each object in Tables [8] and [9] Estimates of 

P[rocky] using an, alterato density model proposed re- 
cently by [Rogers lead to similar assessments, at 
least e come 

'The results of our calculations indicate that all of the 
planets in our sample have a ~50% or greater chance 
of being in the temperate region of their parent stars. 
KOI-0571.05, 1422.04, and 4742.01 are the best candi- 
dates for residing within the HZ boundaries, with P[HZ] 
equal to 98.496, 99.796, and 96.996, respectively. The first 
two receive roughly one third of the flux of the Earth, 
while the latter receives two thirds. Based on their mea- 
sured sizes, KOI-0571.05 and KOI-4742.01 have the best 
chance of being rocky; KOI-1422.04 is larger and has a 
somewhat smaller chance. KOI-3284.01 is the object in 
our sample with the highest probability of being rocky 


in nature (69.6%); its size is consistent with being the 
same as that of the Earth: R, = 1.12*019 Re. Other 
planets with non-negligible chances of being rocky are 
KOI-1422.05, 4087.01, 4622.01, and 4742.01 On the 
other hand, the larger planets KOI-2529.02, 4005.01, and 
4745.01 have [rocky] ~ 5% or smaller, and are most 
likely mini-Neptunes. KOI-3255.01 may still be rocky 
despite its radius of 2.14 Rg, with [rocky] = 11.7%. 

In terms of their insolation, four of our validated plan- 
ets with a sizable chance of being rocky have Sg values 
within 5096 of that of the Earth: KOI-1422.05, 3284.01, 
4087.01, and 4742.01. Of these worlds, KOI-3284.01 and 
KOI-4742.01 appear to be the most “Earth-like” when 
considering both their size and amount of incident flux 
jointly, surpassing KOI-0571.05, which is also rocky but 
receives only one third of the incident flux of our planet. 


30 KOI-4427.01 would otherwise join this group, but it is not 
formally validated as a planet to the same level as the others. Its 
insolation is only about a quarter of that of the Earth. 
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10. DISCUSSION AND CONCLUDING REMARKS 


Eleven of the twelve candidates in our sample have 
been validated here as planets to a very high degree of 
confidence exceeding 99.73% (3e). The high bar we have 
set in this work is commensurate with the importance 
we assign to these planets, which are all small and po- 
tentially in the HZ of their host stars. Few such objects 
have been verified so far, and our work now essentially 
doubles the number of cases listed in Sect. [I] 

'Three of our targets, KOI-0571.05, KOI-1422.04, and 
KOI-1422.05 are in systems of five planets each, and have 
been previously validated by others using very different 
methodologies and assumptions. Our BLENDER valida- 
tions confirm those assessments using conservative and 
realistic hypotheses, generally with a higher degree of 
confidence. KOI-1422.04 and KOI-1422.05 were included 
in the studies of [Lissauer et al] and 
(2014), who considered them to be validated based on 
the general conclusion that the vast majority (> 99%) 
of Kepler candidates in multiple systems are likely to 
be bona-fide planets, as they showed statistically] In 
the present work we are able to assign a specific con- 
fidence level (99.74% and 99.89%) to these two plan- 
ets using BLENDER. KOI-0571.05 was validated separately 
by Oan ana etal] through statistical arguments 
somewhat similar to those in BLENDER in some respects, 
but not in detail. For example, (2014 
included short period planets in their estimate of the 
planet prior, whereas we chose to count only planets with 
similar period as the candidate, because the rates of oc- 
currence of transiting planets are a strong function of 
orbital period. For a planet with a long period of ~130 
days such as KOI-0571.05, including short periods over- 
estimates PL by about a factor of 50 according to our 
calculations. On the other hand, including short periods 
also for the false positives (both for planets and EBs), 
> live done, vll fend qoare 
the blend frequencies higher. This partly offsets the bias 
from the larger planet prior, but still results in a final 
odds ratio a factor of 2 greater than if short periods are 
not included. An additional difference is that in order 
to reach a comfortable validation level 
applied a “multiplicity boost” to their planet prior 
(increasing it by a factor of 30), invoking the fact that 
KOI-0571.05 is in a five-planet system, that the other 
four planets were considered validated by [Lissauer et al] 
and (2014), and that false positives 
are significantly less common in multis. While this is 
a perfectly legitimate argument, the odds ratio achieved 
with BLENDER is high enough (for all our multis) that this 
boost was not required here. Had we used it, our confi- 
dence level for KOI-0571.05 would be 30 times higher. 

One of our candidates, KOI-4427.01, does not quite 
reach our threshold for validation, although in other 
contexts a 99.1696 confidence level such as we achieve 
might be considered high enough to declare it a bona- 
fide planet. Interestingly, this is the same candidate for 
which our centroid motion analysis indicated less confi- 
dence than in other cases that the transit signal comes 


31 We note that another of our targets, KOI-2529.02, is a candi- 
date in a two-planet system that was also included in the statistical 
studies above, but it was not considered to be validated there due 
to unresolved issues with the follow-up observations. 


from the target itself as opposed to a nearby background 
location, consistent with indications from BLENDER. We 
note also that the host star does not show a clear ro- 
tational signature, which causes its properties to be less 
well determined than others for which we were able to 
use an age estimate from gyrochronology to strengthen 
the model fits. For these reasons we do not include this 
object among our formally validated planets, although it 
is entirely possible that additional follow-up observations 
may enable a higher level of validation in the future. It 
remains a very interesting object, as it is likely to be 
within the HZ and also has a fair chance of being rocky. 

We find that the mean densities derived photomet- 
rically from the individual planets in each of the five- 
planet systems KOI-0571 and KOI-1422 show very good 
agreement, supporting the notion (for the first time in 
KOI-1422) that in each case the five planets orbit the 
same star. This is further supported by a similar con- 
sistency found among the limb-darkening coefficients (u1 
and u3), which were fitted for separately for each of the 
planets. For the four targets in our sample with close 
companions the location of the planets is more uncer- 
tain, however. In one case (KO1-2529.02) both BLENDER 
and AP strongly favor the planet transiting the primary 
star. For KOI-1422.04, KOI-1422.05, and KOI-3284.01 
the evidence is in the same direction but not as strong, 
and for KOI-3255.01 both BLENDER and AP are inconclu- 
sive as to which star has the planet. In these four cases 
even if the planet orbits the companion the inferred radii 
are still relatively small (1.68-2.25 Rẹ), and at least for 
KOI-1422.05 and KOI-3284.01 there is still a good chance 
they are rocky. 

Except for KO1-2529.02, the inferred radii of all our 
validated planets are nominally smaller than 2.4 Rg. In 
some cases they are rather different than those reported 
previously, either because of differences in the stellar 
radii, or because some of the previous studies have not 
accounted for dilution by nearby companions. The stellar 
properties in this work are all based on spectroscopic esti- 
mates of the temperature and metallicity, as well as log g 
for the hotter stars and other constraints on luminos- 
ity for the cooler ones, superseding properties reported 
previously that have relied only on photometry. 

'This study has significantly increased the number of 
small validated transiting planets (< 2.5 Rẹ) that are 
potentially in the HZ of their host star. Excluding 
those in our sample that have been announced pre- 
viously, as well as KOI-4427.01, we now add to the 
list KOI-1422.05(*), KOI-2529.02, KOI-3255.01(*), KOI- 
3284.01(*), KOI-4005.01, KOI-4087.01(*), KOI-4622(*), 
KOI-4742.01(*), and KOI-4745.01. The ones followed by 
an asterisk are small enough to be rocky, with a some- 
what reduced chance of that being the case for KOI- 
3255.01. We note that rather than relying on an esti- 
mate of the equilibrium temperature for our assessments 
regarding the HZ, here we have used the effective insola- 
tion, which, as pointed out by [Kopparapu et al] (2013), 
dispenses with the uncertainty in the albedo. Our HZ 
statements are quantified by providing a probability that 
the planet lies within the adopted inner and outer edges 
of the region, which properly accounts for the uncer- 
tainties in all measured stellar and planetary properties. 
Similarly, we have provided a quantitative indication of 
the probability that each of these planets is rocky, also 
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accounting for uncertainties. 

A number of the previously announced small HZ plan- 
ets (see Sect.[I) have radii in excess of 2 Rg, and accord- 
ing to our prescription for P [rocky] they are significantly 
less likely to be rocky than perhaps deemed to be the case 
in the original publications (particularly Kepler-22 b and 
Kepler-296f, with radii of 2.38 Rg and 2.31 Rg, respec- 
tively). Thus, with the new examples added here, the list 
of rocky planets in the HZ is enlarged by a factor of two. 
KOI-3284.01 and KOI-4742.01 are now the validated, 
rocky, HZ planets that appear most similar to the Earth 
when considering both their size (Rp = 1.12*0:19 Re and 


Ry = 1.34*011 Re) and insolation (Ses = 1.40*0 97 Sq 
and Seg = 0.667071 Sg) jointly. 
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